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ABSTRACT 


Formulation  and  application  of  a  three-dimensional  compressible 
turbulent  boundary-layer  analysis  is  presented  for  subsonic  and  tran¬ 
sonic  flow  over  a  yawed  airfoil  of  infinite  extent.  The  governing  tur¬ 
bulent  boundary- layer  equations  are  integrated  using  an  implicit  finite- 
difference  procedure  in  conjunction  with  a  scalar  eddy  viscosity  model 
of  three-dimensional  turbulence.  Comparisons  with  other  analysis 
techniques  as  well  as  experimental  measurements  under  subsonic  wind 
tunnel  conditions  are  presented  to  establish  and  ascertain  the  basic 
validity  and  applicability  of  the  current  technique.  Also  considered 
are  the  effects  of  a  hot  wall  on  the  transonic,  three-dimensional,  tur¬ 
bulent  boundary  layer  which  have  practical  application  to  transonic 
Space  Shuttle  reentry,  where  the  wing  surface  temperature  may  reach 
soak  values  on  the  order  of  twice  the  free- stream  stagnation  tempera¬ 
ture  because  of  the  hypersonic  high- heating  phase  of  the  reentry  tra¬ 
jectory.  Ground  testing  of  Space  Shuttle  configurations  under  contin¬ 
uous  transonic  flow  conditions  with  an  adiabatic  wall  may  not  be  totally 
applicable  to  actual  Shuttle  entry;  airfoil  static  stall  angle,  maximum 
lift  coefficient,  local  skin-friction,  and  other  turbulent  boundary- layer 
parameters  are  strongly  affected  by  hot- wall  relative  to  adiabatic- wall 
conditions. 
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SECTION  I 
INTRODUCTION 


The  three-dimensional  compressible  turbulent  boundary  layer  on 
an  infinite  yawed  airfoil  under  subsonic  or  transonic  flow  conditions  is 
a  problem  which  has  received  little  attention  in  the  literature.  Only 
within  recent  years,  with  the  development  of  accurate  numerical  in¬ 
tegration  methods  using  a  digital  computer,  has  even  the  simplier 
case  of  laminar  compressible  boundary- layer  flow  over  an  infinite 
yawed  wing  become  amenable.  (See  Ref.  1  for  documentation  of  one 
such  analysis. )  Needless  to  say,  the  infinite  yawed  airfoil  never  oc¬ 
curs  in  practical  applications  (e.g. ,  a  transport  aircraft).  The  im¬ 
portance  of  this  type  flow  concerns  the  trends  which  it  reveals  (i.e., 
the  effects  of  three-dimensional  flow  on  boundary- layer  characteristics). 
In  addition,  high  aspect  ratio  finite- length  swept  airfoils  may  be  ade¬ 
quately  represented  by  an  infinite  yawed  wing  in  regions  removed  from 
the  wing  root  and  tip;  an  example  of  such  an  airfoil  is  given  in  Ref.  2 
relative  to  helicopter  rotor  applications. 

The  boundary  layer,  be  it  laminar  or  turbulent,  on  an  infinite 
yawed  wing  is  three-dimensional  in  the  sense  that  the  velocity  and 
shear- stress  profiles  are  skewed,  but  quasi- two-dimensional  insofar 
as  spanwise  derivatives  of  the  pressure,  velocity,  and  temperature  are 
identically  zero.  For  example,  there  is  spanwise  flow  on  the  infinite 
yawed  wing  but  no  net  boundary- layer  "drift"  towards  the  wing  tip  such 
as  occurs  on  finite  airfoils.  An  analytical  study  of  the  complete  three- 
dimensional  turbulent  boundary- layer  flow  over  a  finite  swept  airfoil  is 
an  extremely  difficult  problem,  with  Ref.  3  being  the  only  such  analysis 
currently  available  (to  the  author's  knowledge).  With  respect  to  infinite 
yawed  wing  flows,  many  previous  investigations  (Refs.  4  through  10) 
have  developed  analytical  techniques  applicable  under  incompressible 
conditions  for  both  laminar  and  turbulent  boundary  layers.  For  the 
case  of  compressible  boundary- layer  flow  over  an  infinite  yawed  air¬ 
foil,  only  the  studies  of  Ref.  1  for  laminar  flow  and  Ref.  11  for  turbu¬ 
lent  flow  have  been  performed. 

The  present  report  will  examine  three-dimensional  compressible 
turbulent  boundary- layer  flow  over  an  infinite  yawed  airfoil  under  sub¬ 
sonic  and  transonic  conditions.  The  governing  turbulent  boundary- 
layer  equations  are  integrated  using  an  implicit  finite- difference  pro¬ 
cedure  in  conjunction  with  a  scalar  eddy  viscosity  model  of  three- 
dimensional  turbulence.  Comparisons  of  the  present  approach  with 
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experimental  data  are  presented  in  order  to  establish  and  ascertain  the 
basic  validity  and  applicability  of  the  current  technique  under  subsonic 
and  transonic  flow  conditions. 


SECTION  II 

ANALYTICAL  ANALYSIS 


2.1  OVERVIEW 

The  present  analytical  investigation  employs  a  three-dimensional 
turbulent  boundary- layer  analysis  based  on  a  scalar  eddy  viscosity 
model  of  turbulence  for  the  case  of  subsonic  and  transonic  flow  over 
an  infinite  extent  yawed  airfoil.  Full  development  of  the  boundary- 
layer  analysis  will  be  presented  below;  the  inviscid  analysis  utilizes 
existing  well -documented  procedures  which  are  described  briefly. 


2.2  GOVERNING  BOUNDARY-LAYER  EQUATIONS 

The  present  analysis  employs  the  three-dimensional  compressible 
turbulent  boundary- layer  equations  in  terms  of  time- averaged  mean- 
flow  quantities  as  derived  by  Vaglio-Laurin  (Ref.  12).  Assuming  the 
ratio  of  boundary- layer  thickness  to  local  surface  curvature  to  be  every 
where  small,  the  governing  equations  of  motion,  in  terms  of  the  ortho¬ 
gonal  coordinate  system  x,  y,  z  illustrated  in  Figs.  1  and  2,  Appendix 
I,  reduce  to  the  following  (see  NOMENCLATURE  for  terminology): 


CONTINUITY 


0 


x- MOMENTUM 

O 

_ diT  - <5 Jr 

* f,v r, +  '’“it  - 

y- MOMENTUM 


d  f  <9ii 
dy  If  dy 


(1) 

(2) 

(3) 
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z- MOMENTUM 


ENERGY 


o  o 

—  <9w  —  dyf  d  f  - 7—7 1 

P*d7+  Pwd v  +  pwif  =  -&  +  Ty)Pd7  -  P  VWJ 

o 

—  as  as  _ 3/  a  r  ('as  1  -  p>  as\  -  -***1 

f'  ar  +  cvaT  *  '"Jr  ‘  57 rw  +  p  ,h  J 


where 


V  =  v  +  *-= 
P 


(4) 


(5) 


(6) 


H  =  F 


—  2  -2 
u  +  w 


(7) 


and  the  usual  expressions  for  the  mean  and  fluctuating  parts  of  the  de¬ 
pendent  variables  are  used;  e.  g. , 

P  **  p  +  p'  (8) 

Implicit  in  Vagi io-Laur in’s  derivation  of  the  above  equations  are  the 
following  stipulations: 

a.  The  rates  of  change  of  the  mean-flow  properties  in  the  x-  and 
z-directions  [0(1)J  are  smaller  than  the  rates  of  change  in  the 
y-direction  [0(6"^)]  by  an  order  of  magnitude. 

b.  Mean  squares  and  products  of  the  turbulent  fluctuations  are 
0(6);  that  is,  the  turbulent  level  is  small.  The  terms  involv¬ 
ing  mean  squares  of  the  velocity  fluctuations  are  taken  to  be 
negligible,  which  is  valid  for  high  Reynolds  number  flows  with 
a  zero  or  favorable  pressure  gradient. 

c.  The  time- average  molecular  transports  are  approximated  by 
those  pertaining  to  the  mean  flow  properties;  indeed,  even  the 
latter  are  negligible,  except  very  near  the  wall,  compared  with 
terms  involving  the  turbulent  transports. 

Also  implicit  in  the  above  equations  is  the  requirement  of  an  infinite  ex¬ 
tent  body  of  the  yawed  cylinder  or  yawed  wing  type,  which  leads  to  the 
term  3  /3z  =  0  as  discussed  in  Section  I. 
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If  subscript  w  denotes  wall  and  subscript  e  denotes  outer  .edge  of 
the  boundary  layer,  the  associated  boundary  conditions  on  the  above 
defined  equations  are 

MOMENTUM 

:u  =  v=  w=  uv  -  vw  =  p  v  =  0 


as  y  -»  oo  :  u  -*  U  ,  w  -»  ff . 

e  (9) 

u'v'  -*0,  v'w'  -*  0  ,  p 'v '  -*  0 

ENERGY 

y  «=  0:  H  =  Hw  =  hw,  vU'x  0  (Prescribed  Wall  Enthalpy) 

y  =  0:  -3^-  =^1  =  0,  v'H '  =  0  (Adiabatic  Wall) 
ay  dy 


as  y  -»  oo:  H-*He,  v'H*  -»  0 


(10) 


which  reflect  the  requirements  of  no  slip  and  no  mass  transfer  (suction 
or  blowing)  at  the  wall,  as  well  as  a  prescribed  wall  enthalpy  or  an  ad¬ 
iabatic  wall.  The  y-momentum  equation  (3)  reveals  that  the  static  pres¬ 
sure  variation  across  the  boundary  layer  is  negligible,  and  hence  the 
static  pressure,  p(x),  is  regarded  as  an  external  input  to  the  boundary- 
layer  analysis  from  a  separate  inviscid  analysis.  The  outer  edge  veloc¬ 
ities,  Ue  and  We,  as  well  as  the  outer- edge  static  enthalpy,  he,  must 
be  determined  from  the  inviscid  analysis  consistent  with  the  imposed 
static  pressure  distribution. 

The  gas  model  adopted  for  the  present  study  is  thermally  and  cal- 
orically  perfect  air  having  a  constant  specific  heat  ratio  7  =  1.  40  and 
obeying  the  equation  of  state 


p  o  pRT  (11) 

where  R  =  1716  ft^/sec^-°R.  Under  this  assumption  the  static  enthalpy, 
h,  is  given  by 

h  =  CpT  (12) 

where  Cp  =  6006  ft^/sec^-°R.  The  laminar  viscosity,  p,  is  taken  to 
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obey  Sutherland’s  law 


2.270  x  10~8  T3/2  Ibl-sec 
T  +  198.6  ft2 


(13) 


where  T  must  have  units  of  °R.  The  laminar  Prandtl  number,  Pr,  is 
taken  to  be  a  constant  value  of  0.  71  across  the  entire  boundary  layer: 
i.e.. 


Pr  =  0.71 


(14) 


2.3  TURBULENT  TRANSPORT  MODEL 

Before  Eqs.  (1),  (2),  (4),  and  (5)  can  be  solved,  expressions  must 
be  supplied  for  the  Reynolds  stress  or  turbulent  shear  terms  in  the  mo¬ 
mentum  equations  and  the  turbulent  flux  of  total  enthalpy  in  the  energy 
equation.  The  approach  used  in  the  present  analysis  is  to  model  these 
terms  as  functions  of  the  mean- flow  variables  following  Adams  (Refs. 
13,  14,  and  15),  whose  studies  are  based  on  the  original  analysis  by 
Hunt,  Bushnell,  and  Beckwith  (Ref.  11). 

The  concept  that  the  Reynolds  stress  in  turbulent  flow  is  propor¬ 
tional  to  a  momentum  exchange  coefficient  times  the  mean- flow  velocity 
gradient  normal  to  the  surface  is  well  known  and  commonly  used  in  tur¬ 
bulent  boundary- layer  analyses.  This  concept  is  based  on  an  assumed 
analogy  between  the  so-called  eddy  viscosity  and  the  molecular  viscos¬ 
ity.  The  total  shear  components  in  the  x-  and  z- directions  are  written 
as 


c?u 

_ _ 

du 

du 

»Ty 

-  p  U  '  V  ' 

=  "  Ty 

+  f*  Ji 

(15) 

d  w 

—  *  * 

dw 

(Jw 

11  dy 

—  p  V  w 

"  *37 

+  d7 

(16) 

where  the  eddy  viscosities  ex  and  ez  might,  in  general,  be  different. 
Applying  the  Prandtl  mixing- length  hypothesis  in  conjunction  with  the 
assumption  that  the  eddy  viscosity  is  ,a  scalar  function  independent  of 
coordinate  direction  (which  means  physically  that  the  turbulent  shear 
stress  acts  in  the  mean  rate  of  strain  direction)  results  in  an  eddy  vis¬ 
cosity  relationship  of  the  form 


f 


(17) 
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where  G  is  a  scalar  velocity  function  defined  by 


dC 

dy 


(18) 


The  quantity  i*  is  termed  the  mxing  length  and  is  some  characteristic 
length  related  to  the  size  or  scales  of  eddies  responsible  for  the  flux  of 
momentum  in  the  y- direction.  Under  the  above  model  the  turbulent 
shear  stress  in  a  three-dimensional  turbulent  boundary  layer  may  be 
written  in  the  form 


y  *  Gt)2]  “ 

The  complete  derivation  of  Eqs.  (17),  (18),  and  (19)  is  given  in  Refs. 
13,  14,  and  15. 


The  expression  for  the  total  heat  flux  in  a  turbulent  boundary  layer 
may  be  written  in  terms  of  the  static  enthalpy  as 


JL 

CP  dy 


p  v'h ' 


k  dh 


K  dh 

CP  dy 


(20) 


where  k  is  the  laminar  (molecular)  thermal  conductivity  and  *  is  the 
so-called  eddy  thermal  conductivity.  Using  the  definition  of  the  laminar 
(molecular)  Prandtl  number 


Pr  = 


(21) 


and  defining,  by  analogy,  a  turbulent  Prandtl  number  (based  on  the  use 
of  static  enthalpy)  as 


P't  = 


(22) 


with  e  the  eddy  viscosity  discussed  previously,  the  total  heat  flux  ex¬ 
pression,  Eq.  (20),  may  be  written  in  the  form 


q  = 


f  1  dji 
+  p  Pr(  dy 


(23) 
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2.4  MIXING-LENGTH  MODEL 


The  turbulent  shear  stress  in  a  three-dimensional  turbulent  bound¬ 
ary  layer  as  governed  by  Eq.  (19)  is  treated  herein  by  the  use  of  a  two- 
layer  inner-outer  model  using  Prandtl's  mixing- length  hypothesis  and 
a  modification  of  van  Driest' s  analysis  for  the  near- wall  region.  This 
results  in  a  continuous  distribution  of  the  shear  stress  from  the  lami¬ 
nar  value  at  the  wall,  through  the  fully  turbulent  region,  reaching  zero 
at  the  outer  edge  of  the  boundary  layer.  The  energy  transport  in  a 
turbulent  boundary  layer  is  treated  in  this  work  through  the  incorpora¬ 
tion  of  the  eddy  conductivity,  k,  into  the  turbulent  Prandtl  number,  Pr^.. 

In  the  manner  of  Escudier  (Ref.  16),  Patankar  and  Spalding  (Ref. 
17)  recommend  the  following  variation  of  the  mixing  length,  SL+,  across 
the  turbulent  two-dimensional  boundary  layer  which  is  adopted  for  the 
present  three-dimensional  case  by  noting  that  the  scalar  properties  of 
a  turbulence  field  are  unlikely  to  be  affected  by  moderate  three-dimen¬ 
sionality  because  turbulence  is  inherently  three-dimensional  in  nature 
for  even  so-called  two-dimensional  flows: 

£.  =  k*y,  for  0  <  y  <  Xyj/k* 

(24) 

K  -  Ayfi  h>r  *yp/k,  <  y 


where  the  values  for  the  various  numerical  constants  are  taken  to  be 
k*  =  0.435  and  X  =  0.09.  The  value  of  y  at  the  point  where  the  velocity 
in  the  boundary  layer  is  equal  to  0.99  of  the  velocity  at  the  boundary- 
layer  outer  edge  is  used  to  define  the  distance  y^;  i.e.. 


y  £ 


2  2-1 

,  ,  [ro  +  (w)i 

y- value  where  ■ 

*  1/2 
[(ue)2  +  <we)2] 


0.99 


(25) 


By  analogy  with  Stokes'  solution  for  an  infinite  flat  plate  undergoing 
simple  harmonic  motion  parallel  to  itself  in  an  infinite  fluid,  van  Driest 
(Ref.  18)  concluded  that  in  the  vicinity  of  a  wall  the  total  shear  stress  in 
a  turbulent  two-dimensional  fluid  should  be  of  the  form 


(26) 
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which  results  in  an  exponential  damping  of  the  turbulent  part  of  the 
shear  stress  as  the  wall  is  approached  and  yields  exactly  the  laminar 
shear  stress  form,  t=  iu( 9u/3y),  at  the  wall.  Although  Eq.  (26)  was 
originally  developed  for  incompressible  flow,  it  can  be  applied  to  com¬ 
pressible  flow  by  application  of  the  suggestion  by  Patankar  and  Spalding 
(Ref.  17)  that  the  local  value  of  shear  stress  be  used  instead  of  the  wall 
value  as  originally  recommended  by  van  Driest  (Ref.  18).  Hence,  by 
analogy  of  Eq.  (26)  with  Eqs.  (18)  and  (19),  the  relationship  for  the 
three-dimensional  near- wall  shear  stress  as  used  in  the  present  anal¬ 
ysis  is 


T 


(27) 


where  the  constant  A*  is  taken  to  be  26.  0  following  the  original  van 
Driest  proposal  (Ref.  18).  Note  that  the  damping  term  in  Eq.  (27)  re¬ 
flects  the  application  of  the  local  total  shear  stress  as  opposed  to  the 
wall  shear  stress  of  Eq.  (26),  as  discussed  previously. 

Based  on  Eqs.  (17),  (18),  (24),  and  (27),  the  eddy  viscosity  expres¬ 
sion  for  the  inner  region  is 


dy 


(28) 


and  for  the  outer  region  is 


( 


O 


p\z 


2  dG 

ne7 


(29) 


with  the  constants  k#.  A#,  X,  and  y^  defined  previously.  The  constraint 
used  to  define  the  end  of  the  inner  region  and  the  beginning  of  the  outer 
region  is  the  continuity  of  the  eddy  viscosity.  From  the  wall  outward, 
the  expression  for  the  inner  eddy  viscosity  applies  until  ei  =  eQ,  from 
which  point  the  outer  eddy  viscosity  is  used.  A  schematic  of  this  vari¬ 
ation  in  terms  of  the  mixing  length  follows. 
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Constant 


The  turbulent  Prandtl  number  (based  on  the  static  enthalpy  defini¬ 
tion  of  the  turbulent  heat  flux)  as  given  by  Eq.  (22)  is  physically  a  meas¬ 
ure  of  the  ratio  of  the  turbulent  transport  of  momentum  to  the  turbulent 
transport  of  heat.  For  the  present  work,  the  turbulent  Prandtl  number 
defined  by  Eq.  (22)  is  taken  to  remain  constant  at  the  value  0.  90  across 
the  entire  boundary  layer  as  recommended  by  Patankar  and  Spalding 
(Ref.  17)  for  two-dimensional  turbulent  boundary  layers. 

For  the  case  of  purely  laminar  flow  in  the  present  analysis,  the 
inner  and  outer  eddy  viscosities  defined  by  Eqs.  (28)  and  (29)  are  both 
set  identically  equal  to  zero  so  that  the  governing  turbulent  boundary- 
layer  equations  reduce  to  their  laminar  counterpart. 
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2.5  ALTERNATE  FORM  OF  ENERGY  EQUATION 

To  apply  the  eddy  viscosity- eddy  thermal  conductivity  model  dis¬ 
cussed  above  to  the  present  work,  the  energy  equation  (5)  must  be  re¬ 
written  to  replace  the  total  turbulent  energy  flux  p  v'H'  in  terms  of 
static  quantities.  By  following  the  steps  outlined  on  pages  176-183  in 
Dorrance  (Ref.  19),  one  can  see  that  the  energy  equation  (5)  may  be 
written  in  the  equivalent  form 


O 


based  on  the  definitions  of  the  scalar  eddy  viscosity  and  the  turbulent 
Prandtl  number  given  in  Section  2.  3  of  the  present  report.  Equation 
(30)  is  to  be  regarded  as  the  three-dimensional  analog  of  Eq.  (7-31)  on 
page  182  in  Dorrance  (Ref.  19)  for  a  two-dimensional  turbulent  boundary 
layer.  Implicit  in  the  derivation  of  Eq.  (30)  is  the  neglect  of  terms  in¬ 
volving  mean  squares  of  the  velocity  fluctuations,  which  is  consistent 
with  Vaglio-Laurin's  derivation  of  the  governing  equations  as  discussed 
previously  in  Section  2.  2. 


2.6  COORDINATE  TRANSFORMATION 

To  facilitate  numerical  integration  of  the  governing  boundary- layer 
equations  (1),  (2),  (4),  and  (30),  it  is  convenient  to  transform  them  to 
a  coordinate  system  that  removes  the  mathematical  singularity  at  x  =  0 
and  stretches  the  coordinate  normal  to  the  surface,  as  is  usually  done 
in  two-dimensional  laminar  flow  analyses.  The  coordinate  transforma¬ 
tion  used  in  the  present  work  is  the  well-known  Illingworth- Levy  (Ref. 
19,  p.  30)  transformation.  The  new  independent  variables  introduced 
are 


f(x)  =  /  pe/ieUe dx 

0 

(31) 

rj(x,y)  =  f  j-  dy 

\/  o  Pt 

(32) 
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so  that  the  transformed  streamwi.se  and  normal  derivatives  become 


d  tT  d  drt  d 

Tx  =  5?  +  dx  d^i 

d  d 

dy  =  vW~  drf 


(33) 

(34) 


Define  a  so-called  stream  function,  p(x,y),  in  such  a  manner  as  to  iden¬ 
tically  satisfy  the  continuity  equation  (1);  i.e.. 


pu  = 


dp 

dy 


(35) 


p\  = 


dp 

dx 


(36) 


and  introduce  a  nondimens ional  stream  function  f(£,  rj)  such  that 

p(^,r})  =  \/W  f  (£.»?)  (37) 

so  that  the  governing  boundary- layer  equations  (2),  (4),  and  (30)  become, 
in  the  transformed  (£,  rj)  coordinates: 

x- MOMENTUM 


<?Yr  +  ff"+  &6-(r)2]  =  -  f"^] 


z- MOMENTUM 


ENERGY 


(fcT  *  fc- .  2f[r§|  -  c- *] 


(£»)'♦  f«'  +  O'  - 


with  the  new  dependent  variables 


*'<£,)  =  f- 

u  e 


(38) 


(39) 


(40) 


(41) 
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c<£>?)  -  ^  (42) 

8  -  ip  (43) 

e 


The  following  definitions  apply  to  the  above  equations: 


Pe  T 

p  = 

(44) 

It- 

PePe 

(45) 

P 

(.46) 

(47) 

it  dU<= 
ue“dT 

(48) 

(49) 

where  the  inviscid  x- momentum  equation  evaluated  on  the  body  surface. 


dp 

T*  +  P 


d  0. 


e  Ue 


dx 


(50) 


has  been  used  to  relate  the  inviscid  velocity,  Ue,  to  the  imposed  invis¬ 
cid  static  pressure,  p.  In  the  transformed  governing  equations,  primes 
denote  partial  differentiation  with  respect  to  the  rr  coordinate;  i.e.. 


f' 


r 


(51) 


An  expression  for  the  static  temperature  ratio,  9,  of  Eq.  (44)  can  be 
found  by  application  of  Eq.  (7)  to  yield 
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0  “  T"  ■  r 

e  e 


5  -  (4^) 


(52) 


which,  using  Eqs.  (41),  (42),  and  (43),  can  be  written  in  the  form 

2  «,  2 


Q 


.  g  _  I  \—  (f')2  +  —  (c)2 
he  B  2  |_he  he  KC> 


(53) 


The  physical  boundary  conditions  given  by  Eqs.  (9)  and  (10)  become, 
in  terms  of  the  transformed  variables, 


MOMENTUM 


ftf.Tj  -  0)  =  0 

{‘(fa  -  0)  «  0 

c(£jj  -  0)  -  0 

iim  f'(f,j j)  -  1 

lim  c  (£,>?)  -  1 

7f-*co 


(54) 


ENERGY 


«w  hw 


g(£. v  =  0)  =  jj— ■  =  jj—  =  gw  (Prescribed  Wall  Enthalpy) 
g'(£»?=  0)  =  0  (Adiabatic  Wall) 


=  1- 


(55) 


2.7  BOUNDARY-LAYER  PARAMETERS 

Given  the  numerical  solution  to  the  governing  equations  of  motion 
(Eqs.  (38),  (39),  and  (40))  following  the  integration  procedure  of  Appen¬ 
dix  III,  the  associated  local  boundary-layer  parameters  at  a  given  body 
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station  may  be  determined  as  follows.  The  local  convective  heat  flux 
at  the  body  surface  (y  =  0)  is  given  by  the  well-known  Fourier  law, 


which  can  be  written  in  the  equivalent  form 


w 


p*(d h\ 

Pr  \<3y  / 


making  use  of  Eqs.  (7)  and  (21)  with 

/dH\  / dh\  \Jfc 

U/w  “  Ww  +\f  d y  V  <Jy_|w 


(56) 


(57) 


(58) 


through  application  of  u  =  0  and  w  =  0  at  the  wall,  as  given  by  Eq.  (9). 
In  terms  of  the  transformed  (§,  rj)  coordinates,  Eq.  (57)  becomes 


where 


—  £  p  u  u  H 

qw  - - — -  g 

pWW 


0) 


t 


w 


P  w^w 
P'Pe 


(59) 


(60) 


In  a  similar  manner,  the  local  wall  shearing  stress  in  the  x-  and 
z-directions  may  be  written  in  terms  of  the  transformed  coordinates  as 


r 


W,  X 


r 


w,z 


y/W 


£w  P.f*.  U.W.  „  t 

- — -  C  (£lj 

y/W 


0) 


0) 


(61) 

(62) 


Physical  height  in  the  boundary  layer  is  determined  from  the  y- 
transformation  relationship,  Eq.  (34),  as 
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V-jFf-  /’e*l  (63) 

r  e  u  e  O 

For  the  case  of  three-dimensional  boundary- layer  flow,  two  lengths 
characterizing  mass-flow  defect  may  be  defined  in  terms  of  the  profiles 
of  the  two  velocity  components  in  the  x-  and  z- directions 


-  f  (‘  - 


In  terms  of  the  transformed  (?,  rj)  coordinates,  the  above  two  mass-flow 
defect  lengths  may  be  written  as 


"e  e  o 


As  shown  by  Moore  (Ref.  20,  p.  7),  Eq.  (64)  is  the  proper  definition  of 
the  boundary-layer  displacement  thickness  for  flow  over  an  infinite  ex¬ 
tent  yawed  body  (airfoil  or  cylinder);  i.e.,  6*  defined  by  Eq.  (64)  prop¬ 
erly  describes  the  extent  to  which  the  nonviscous  flow  is  deflected  by 
the  boundary  layer  on  an  infinite  yawed  airfoil. 

In  a  similar  manner,  the  momentum  thicknesses  in  the  x-  and  z- 
directions  may  be  written  as 


■  Cik  k  - 1> 


iH_  r~  f'(i-r)  d7 

Pe  0 
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and 


e 


m ,  z 


.  yzT  f“ 

Pe  Ue  o 


c(  1  -  c)  di? 


(70) 


Along  the  stagnation  or  attachment  line  of  a  blunt-nosed  infinite  ex¬ 
tent  yawed  airfoil  where  ?  =  0  since  x  =  0  (see  Fig.  1  for  clarification), 
the  above  equations  reduce  to  the  following  limiting  forms 


(under  the  restriction  that  Ue 


x  near  x  =  0  so  that 


-  1  F d  uel  2 

2  p^ehrjx=0*  nearx  =  0): 


(71) 


(72) 


(73) 


(74) 


/  (0  -  c)  irf 


x=  o 


(75) 


W..J 


3£=  0 


/“f'd-n  a n 


(76) 
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KJ 


Pe 


/  c(l  -  c)  diy 


(77) 


Furthermore,  the  pressure  gradient  parameter,  0,  defined  by  Eq.  (48) 
assumes  the  value  )3  si  on  the  stagnation  or  attachment  line  under  the 
above  restriction  on  Ue  and  5  near  x  =  0. 

Various  quantities  of  interest  are  defined  as  follows: 

FREE- STREAM  REYNOLDS  NUMBER  BASED  ON  CHORD  LENGTH  C 

P«v«c 


Re 


oo,C 


(78) 


LOCAL  UNIT  REYNOLDS  NUMBER  BASED  ON  INVISCID 
EDGE  CONDITIONS 


[Ree/ftl  = 


[Rec/ft], 


PeU« 


Pe  "e 


(79) 

(80) 


LOCAL  MOMENTUM  THICKNESS  REYNOLDS  NUMBER  BASED  ON 
INVISCID  EDGE  CONDITIONS 


Re 


c,8 


p  U  6 
"e  e  m,x 


m,x  Pe 


Re 


e,6_ 


Pe 


Pe 

FREE-STREAM  MACH  NUMBER 

V 


(81) 

(82) 


M. 


vPX 

LOCAL  MACH  NUMBER  BASED  ON  INVISCID  EDGE  CONDITIONS 

U_ 


(83) 


M 


e.x 


*  \/Y  R  Te 


(84) 
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N/y  R  Te 


LOCAL  ANGULAR  DIRECTION  OF  A  STREAMLINE  RELATIVE  TO  THE 
Z-AXIS  DIRECTION 


<u  =  arc  tan 


which  becomes  at  the  outer  edge 


toe  =  arc  tan 


and  at  the  body  surface 


<j3  =  arc  tan 


[Ue  f"  (£17  •  0)“| 


LOCAL  STANTON  NUMBER  BASED  ON  FREE- STREAM  CONDITIONS 


~  Pcc  v«S(To^-Tw> 

LOCAL  SKIN -FRICTION  COEFFICIENT  BASED  ON  FREE- STREAM 
CONDITIONS 


•*  «PooV» 


00  <9 

LOCAL  SKIN- FRICTION  COEFFICIENT  BASED  ON  INVISCID  EDGE 
CONDITIONS 

r  -  f|>ix 
fe.x  'A  pe  uf 


‘  f  =  9 

C,Z  «  Pe  We 


18 


AEDC-TR-73-1 1 2 


LOCAL  BOUNDARY- LAYER  SHAPE  FACTOR 
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f,x 
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f,z 


(94) 

(95) 


2.8  BOUNDARY-LAYER  TRANSITION  CORRELATION  PARAMETERS 


As  discussed  by  Cumpsty  and  Head  (Ref.  6)  the  attachment  line 
flow  on  an  infinite  swept  wing  resembles  fully  developed  pipe  or  channel 
flow  much  more  closely  than  normal  boundary- layer  flow  in  that  the  at¬ 
tachment  line  boundary  layer  is  constant  in  thickness  (and  indeed  in  all 
its  properties)  with  downstream  distance;  see  Fig.  3  for  illustration. 

In  particular,  the  attachment  line  boundary  layer,  be  it  laminar  or  tur¬ 
bulent,  can  be  specified  by  the  value  of  a  single  nondimensional  param¬ 
eter,  C*,  defined  as  (see  Ref.  5  for  clarification) 


C*  = 


(96) 


with  the  nomenclature  defined  as  in  Fig.  4. 


The  experimental  work  of  Cumpsty  and  Head  (Ref.  6)  shows  quite 
clearly  that  the  attachment  line  laminar  boundary  layer  appears  very 
stable  with  respect  to  small  disturbances.  To  produce  a  turbulent  at¬ 
tachment  line  flow,  there  must  be  a  sufficient  disturbance  and  the  pa¬ 
rameter  C*  must  be  large  enough.  Assuming  a  sufficient  disturbance 
(for  example,  a  tripping  device  such  as  a  trip  wire),  the  attachment 
line  boundary  layer  has  the  properties  of  a  fully  turbulent  flow  for  val¬ 
ues  of  C*  greater  than  1.  4  x  10*.  Without  external  disturbance,  the 
laminar  attachment  line  boundary  layer  exhibits  no  instability  charac¬ 
teristics  for  values  of  C*  as  large  as  3.  7  x  10^.  Furthermore,  at  C*  = 
3.  1  x  105  disturbances  decay  from  a  trip  wire  almost  large  enough  to 
produce  immediate  transition  to  a  fully  turbulent  flow.  These  results 
strongly  suggest  that,  at  large  values  of  C*  in  conjunction  with  a  trip¬ 
ping  device  only  slightly  smaller  than  a  "critical"  size,  the  attachment 
line  flow  will  revert  to  a  laminar  condition;  for  a  tripping  device  just 
above  the  critical  size  there  appears  to  be  an  immediate  transition  to 
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fully  turbulent  flow  with  the  downstream  flow  essentially  independent  of 
the  nature  of  the  disturbance.  The  criterion  proposed  by  Gaster  (Ref. 

21)  appears  applicable  and  accurate  for  estimating  the  critical  size  of 
the  tripping  device. 

Assuming  that  the  attachment  line  boundary  layer  has  not  been 
tripped  as  discussed  above,  natural  transition  will  occur  (under  high 
Reynolds  number  conditions)  at  some  downstream  chordwise  location 
because  of  the  three-dimensional  nature  of  the  swept- wing  boundary- 
layer  flow.  Briefly,  three-dimensional  crossflow  has  an  adverse  ef¬ 
fect  on  laminar  boundary-layer  stability  in  that  a  system  of  stream- 
wise  vortices  contained  within  the  boundary  layer  may  be  formed,  ap¬ 
parently  because  of  the  inflection  point  in  the  rotated  crossflow  velocity 
profile  illustrated  in  Fig.  5,  which  is  unstable  to  small  disturbances. 
This  vortex  formation  apparently  signals  the  onset  to  three-dimensional 
crossflow- induced  transition  of  the  boundary  layer  from  laminar  to  tur¬ 
bulent  flow.  A  complete  discussion  of  this- phenomenon  relative  to  sub¬ 
sonic  swept- wing  flows  may  be  found  in  a  report  by  Owen  and  Randall 
(Ref.  22);  supersonic  flows  have  been  examined  by  Chapman  (Ref.  23) 
and  hypersonic  flows  by  Adams  .(Ref.  24).  Application  to  hypersonic 
lifting  bodies  at  incidence  (such  as  the  Space  Shuttle)  have  been  re¬ 
ported  by  Adams  and  Martindale  (Ref.  25). 

The  exact  location  at  which  the  above- discussed  vortex  system  will 
originate  cannot  be  determined  from  classical  boundary- layer  stability 
theory  for  three-dimensional  flows  such  as  that  presented  by  Gregory, 
Stuart,  and  Walker  (Ref.  26).  Instead,  the  abrupt  formation  of  these 
vortices  and  also  the  development  of  complete  turbulence,  i.e.,  transi¬ 
tion,  in  a  three-dimensional  boundary  layer  can  apparently  be  correlated 
with  a  so-called  maximum  local  crossflow  Reynolds  number,  x,  defined 
(Refs.  22,  23,  and  24)  as 


X  = 


Pe  ws£,max  ® 

p7~ 


(97) 


where  w  .  is  the  maximum  crossflow  velocity  in  the  streamline 
sf,max  J 

coordinates  of  Fig.  5,  and  6  is  the  boundary- layer  thickness  defined  as 

the  normal  distance  from  the  surface  where  the  total  resultant  velocity 


—2 
+  W 
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reaches  0.  995  of  the  total  resultant  inviscid  edge  velocity 

vV+W,2  i 

pe  and  jue  are  the  values  of  density  and  viscosity,  respectively,  evaluated 
at  the  inviscid  edge  conditions.  Owen  and  Randall  (Ref.  22)  found  the 
critical  value  of  crossflow  Reynolds  number  for  vortex  formation  and 
for  crossflow- induced  transition  to  be  125  and  175,  respectively,  based 
on  the  leading  edge  of  swept  wings  at  subsonic  speeds.  The  work  by 
Chapman  (Ref.  23)  on  swept  cylinders  at  supersonic  speeds  (free- stream 
Mach  numbers  up  to  seven)  indicates  that 

X<  100 — ►Laminar  Boundary  Layer 

100  ^  X  ^200 — ►Vortex  Formation  arid  Transitional 

Boundary  Layer 

X  >  200— ►Turbulent  Boundary  Layer 

which  means  that  the  crossflow  stability  criterion  of  Owen  and  Randall 
may  apparently  be  expected  to  apply  without  change  on  cylindrical  lead¬ 
ing  edges  for  both  subsonic  and  supersonic  flows.  Chapman's  work  fur¬ 
ther  showed  that  the  amount  of  crossflow  needed  to  induce  crossflow  in¬ 
stability  downstream  of  the  leading  edge  was  very'  small  -  on  the  order 
of  one  to  five  percent  of  the  inviscid  edge  velocity  for  the  conditions 
observed.  This  means  physically  that  transition  may  occur  extremely 
near  the  leading  edge  of  a  swept  blunt-nosed  airfoil  because  of  the  small 
amount  of  crossflow  needed  to  cause  crossflow- induced  transition  at 
small  values  of  the  local  crossflow  Reynolds  number,  X. 

It  should  be  noted  that  the  above-discussed  crossflow- induced  tran¬ 
sition  phenomenon -and  stability  criterion  have  been  verified  under  sub¬ 
sonic  flight  conditions  on  various  aircraft  as  discussed  in  the  Appendix 
of  Ref.  22. 


2.9  NUMERICAL  SOLUTION  OF  THE  GOVERNING  BOUNDARY-LAYER  EQUATIONS 

In  the  present  study  of  infinite  extent  yawed  airfoil  boundary-layer 
flows,  numerical  solution  of  the  governing  nonlinear,  parabolic,  partial 
differential  equations  (Eqs.  (38),  (39),  and  (40))  is  performed  by  ob¬ 
taining  linear  finite-difference  equivalents  of  the  equations  and  solving 
these  using  an  iterative,  marching,  implicit  finite- difference  integration 
technique  involving  inversion  of  tridiagonal  matrices.  Full  details  of 
this  numerical  approach  are  given  in  Appendix  III. 
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The  digital  computer  code  is  written  in  FORTRAN  63  for  use  on  a 
CDC  1604-B  digital  computer.  Solution  time  (including  on-line  input/- 
output)  varies  from  approximately  30  minutes  to  one  hour  for  an  airfoil 
divided  into  100  to  200  stations.  The  implicit  nature  of  the  present  finite- 
difference  approach  has  not  presented  any  numerical  stability  problems. 

With  respect  to  application  of  this  code,  some  operational  details 
for  various  types  of  analyses  are  as  follows.  Given  the  inviscid  flow 
field  over  the  wing  determined  as  in  Section  2.  10,  one  must  choose  the 
grid  mesh  (see  Appendix  III)  and  the  type  of  wall  boundary  condition 
for  the  energy  equation  -  adiabatic  or  prescribed  surface  temperature. 
For  both  subsonic  and  transonic  calculations  under  either  continuous 
wind  tunnel  or  actual  flight  conditions,  the  adiabatic  wall  option  is  nor¬ 
mally  used  since  this  is  the  physically  correct  boundary  condition  for 
this  type  of  flow.  However,  the  surface  temperature  can  be  assigned 
a  constant  value  over  the  entire  wing  if  one  is  interested  in  heat- 
transfer  phenomena;  an  example  of  this  type  analysis  will  be  given  in 
Section  3.  3  relative  to  Space  Shuttle  reentry  applications.  In  the  pres¬ 
ent  report,  unless  otherwise  specifically  stated  (as  in  Section  3.  3),  all 
calculations  have  been  performed  using  an  adiabatic- wall  boundary  con¬ 
dition. 


2.10  INVISCID  FLOW  FIELD 

The  inviscid  flow-field  quantities  needed  for  input  to  the  present 
boundary- layer  analysis  (such  as  Ue,  We,  p,  etc.  )  are  herein  deter¬ 
mined  via  application  of  classical  sweep  theory  for  sweptback  wings 
(Ref.  27,  pp.  415-416).  Briefly,  sweep  theory  applied  to  an  infinite 
yawed  wing  having  sweep  angle  A  (see  Figs.  1  and  2)  is  based  upon  the 
fact  that  the  spanwise  (z-direction)  component  of  the  free- stream  flow, 
namely  V,,,  sin  A,  does  not  influence  the  pressure  distribution  on  the 
wing;  only  the  chordwise  (x- direction)  component,  V0  cos  A,  is  signi¬ 
ficant  for  determination  of  the  inviscid  flow  governing  the  surface  pres¬ 
sure  distribution.  This  means  that  well-established  two-dimensional 
analysis  techniques  can  be  applied  for  determination  of  infinite  yawed 
wing  inviscid  flow  fields  in  the  chordwise  (x)  direction.  For  subsonic 
(incompressible)  flows,  the  surface  pressure  distribution  on  a  specified 
NACA  airfoil  section  can  be  derived  from  the  tabular  data  presented  in 
Appendix  I  of  Abbott  and  von  Doenhoff  (Ref.  28).  For  transonic  (com¬ 
pressible)  flows,  the  surface  pressure  distribution  on  a  specified  air¬ 
foil  section  can  be  determined  from  transonic  small- disturbance  theory 
using  the  digital  computer  codes  of  Refs.  29,  30,  31,  and  32  for 
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numerical  solution  (via  relaxation  procedures)  of  the  transonic  small- 
disturbance  equation  written  in  terms  of  the  velocity  potential.  In  either 
case,  the  airfoil  ordinates  for  NACA  airfoil  sections  can  be  obtained 
from  tabulated  values  in  the  Appendixes  of  Ref.  28.  An  alternate  ap¬ 
proach  for  either  subsonic  or  transonic  conditions  is  an  experimentally 
determined  surface  pressure  distribution;  such  may  require  numerical 
smoothing  before  input  since  numerical  differentiation  of  the  pressure 
distribution  is  used  to  determine  the  local  static  pressure  gradient 
which  is  one  of  the  parameters  required  in  the  boundary- layer  calcu¬ 
lation. 

With  the  surface  pressure  distribution  determined  from  one  of  the 
above  approaches,  all  inviscid  surface  quantities  at  any  chordwise  (x) 
location  along  the  airfoil  can  be  determined  through  application  of  an 
isentropic  expansion  from  the  stagnation  or  attachment  line  conditions 
to  the  desired  x-location  pressure.  Constancy  of  total  energy  is  used 
to  determine  the  x- direction  component  of  velocity,  Ue(x),  knowing  the 
value  of  the  static  temperature,  Te(x),  from  an  isentropic  expansion 
and  the  fact  that  We(x)  remains  constant  at  the  value  We(x)  =  V®  sin  A 
from  sweep  theory.  An  alternate  procedure  useful  for  subsonic  flows 
where  the  inviscid  surface  velocity  in  the  x- direction,  Ue(x),  is  known 
(from  Ref.  28,  for  example)  is  as  follows:  Knowing  the  local  value  of 
Ue(x),  as  well  as  We(x)  =  sin  A  from  sweep  theory,  the  local  value 
of  the  static  temperature,  Te(x),  can.be  determined  from  constancy  of 
total  energy.  The  local  static  pressure  can  then  be  determined  through 
application  of  an  isentropic  expansion  from  the  stagnation  or  attachment 
line  conditions  to  the  local  x- direction  static  temperature.  These  tech¬ 
niques  comprise  a  subroutine  for  determination  of  the  inviscid  edge  con¬ 
ditions  required  in  the  numerical  solution  of  the  boundary- layer  equa¬ 
tions  via  the  digital  computer  code  described  in  Section  2.  9. 


SECTION  III 

RESULTS  AND  DISCUSSION 


Presentation  of  results  from  the  present  analysis  technique  for 
three-dimensional  turbulent  boundary-layer  flow  on  infinite  yawed  bodies 
will  be  given  in  this  section.  Included  are  stagnation  or  attachment  line 
flows  as  well  as  infinite  extent  yawed  airfoils  under  both  subsonic  and 
transonic  flow  conditions.  Comparisons  with  other  analytical  analyses 
and  experimental  measurements  are  presented  in  order  to  ascertain 
the  validity  and  applicability  of  the  current  approach. 
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3.1  STAGNATION  OR  ATTACHMENT  LINE  FLOW 

Certainly  one  of  the  simplest  of  all  infinite  yawed  body  flows  is  the 
classical  yawed  cylinder  case,  which  is  directly  applicable  to  the  stag¬ 
nation  or  attachment  line  flow  on  a  blunt-nosed  infinite  extent  swept  air¬ 
foil.  Cumpsty  and  Head  (Ref.  5)  have  considered  this  problem  in  detail 
for  the  case  of  subsonic  (incompressible)  turbulent  attachment  line  flow 
on  a  blunt-nosed  infinite  extent  swept  wing.  Their  results  can  be  ex¬ 
pressed  in  terms  of  a  single  parameter,  C*,  discussed  previously  in 
Section  2.  8  and  defined  by  Eq.  (96)  as 


with  the  nomenclature  as  in  Fig.  4.  Note  that  C*  may  be  considered 
as  a  stagnation  or  attachment  line  Reynolds  number  based  solely  on 
local  inviscid  conditions;  all  effects  of  sweep,  leading- edge  radius, 
and  free-stream  velocity  are  absorbed  into  the  inviscid  flow  parameters. 

Presented  in  Figs.  6,  7,  and  8  are  comparisons  of  stagnation  or 
attachment  line  calculations. using  the  present  analysis  technique  rela¬ 
tive  to  experimental  measurements  under  subsonic  flow  conditions  as 
reported  by  Cumpsty  and  Head  (Ref.  6).  In  all  three  figures  circles 
denote  experimental  measurements  of  the  spanwise  (z)  direction 
boundary- layer  flow  with  a  0.  064-in.-diam  trip  wire  wrapped  around 
the  leading  edge  about  18  in.  from  the  upstream  tip  of  the  wing;  the  in¬ 
verted  triangles  denote  experimental  measurements  without  the  trip 
wire.  As  discussed  previously  in  Section  2.  8,  the  purpose  of  the  trip 
wire  is  to  provide  a  form  of  upstream  disturbance  necessary  to  produce 
a  turbulent  flow  along  the  attachment  line.  Full  details  of  the  experi¬ 
mental  measurement  techniques  may  be  found  in  Ref.  6. 

The  important  point  to  observe  from  Figs.  6  and  7  is  that  there  is 
no  experimental  evidence  of  any  instability  (i.e.,  tendency  to  undergo 
boundary- layer  transition)  for  the  attachment  line  flow  without  the  trip 
wire.  However,  addition  of  the  trip  wire  causes  an  almost  instanta¬ 
neous  transition  from  laminar  to  fully  turbulent  flow  for  values  of  C* 
greater  than  about  0.6  x  10^  to  0.  8  x  10^.  The  present  numerical  cal¬ 
culations  are  in  reasonable  agreement  with  experiment  for  both  lami¬ 
nar  and  turbulent  flow  conditions;  in  addition,  the  present  calculations 
are  in  good  agreement  with  the  approximate  integral- type  approach  by 
Cumpsty  and  Head  (Ref.  5),  as  can  be  seen  by  comparing  the  results 
in  Figs.  6  and  7  with  Fig.  4  of  Ref.  6.  Further  observe  that  the 
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present  turbulent  calculation  procedure  automatically  reduces  to  a  lami¬ 
nar  flow  calculation  for  small  values  of  the  C*  parameter  (values  of  C* 
less  than  about  0.  2  x  1CP)  with  a  sort  of  calculated  transition  region  be¬ 
tween  C*  values  of  0.  2  x  1(P  and  about  0.  5  x  10^.  This  interesting  be¬ 
havior  can  be  directly  traced  to  the  fact  that  the  present  eddy  viscosity 
model  evaluated  on  the  stagnation  or  attachment  line  is  of  the  form 

M*=o  « 

as  can  be  seen  from  consideration  of  Eqs.  (28)  or  (29)  evaluated  at  x  =  0 
using  Eqs.  (18),  (33),  and  (73)  as  well  as  5  ^  1/2  Peiue(dUe/dx)x=Q  x^ 
near  x  =  0. 

A  comparison  of  the  spanwise  (z-direction)  local  skin-friction  coef¬ 
ficient  on  the  stagnation  or  attachment  line  based  on  the  present  calcu¬ 
lation  technique  is  presented  in  Fig.  8  relative  to  the  experimental 
measurements  (via  Preston  tubes)  reported  in  Fig.  6  of  Cumpsty  and 
Head  (Ref.  6).  The  measurements  are  based  upon  Patel's  Preston  tube 
calibration  using  the  two-dimensional  law  of  the  wall  (see  Ref.  6  for 
details).  In  addition,  the  measurements  correspond  to  a  tripped  up¬ 
stream  flow  such  that  the  boundary  layer  should  be  fully  turbulent  at 
the  measurement  station  location.  In  general,  agreement  of  the  pres¬ 
ent  numerical  calculation  procedure  with  experiment  is  good;  also 
shown  in  Fig.  8  for  the  sake  of  comparison  is  the  so-called  Method  3 
calculation  by  Cumpsty  and  Head  (Ref.  6).  Note  again  the  calculated 
transition  region  (0.  2  x  10^  <  C*  <  0.  5  x  HP)  between  laminar  and  tur¬ 
bulent  flow  as  shown  near  the  top  left-hand  corner  of  Fig.  8. 

All  of  the  above- presented  results  are  for  the  case  of  stagnation 
or  attachment  line  boundary-layer  flow  on  a  blunt-nosed  infinite  extent 
swept  wing  under  low  subsonic  conditions  where  the  flow  is  essentially 
incompressible.  To  the  author's  knowledge  there  has  been  no  work, 
either  analytical  or  experimental,  on  the  transonic  stagnation  or  attach¬ 
ment  line  turbulent  boundary-layer  flow  past  an  infinite  extent  swept 
wing  or  a  yawed  circular  cylinder.  Since  compressibility  effects  be¬ 
come  of  importance  under  transonic  conditions,  it  would  be  of  great 
value  to  conduct  a  joint  analytical /experimental  investigation  of  tran¬ 
sonic  turbulent  boundary  layers  on,  say,  yawed  circular  cylinders  of 
varying  sweep  angle  so  as  to  assess  the  applicability  of  the  C*  similar¬ 
ity  parameter  under  nonsubsonic  flow  conditions. 
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3.2  INFINITE  EXTENT  YAWED  AIRFOIL  UNDER  SUBSONIC  CONDITIONS 

The  case  of  subsonic  incompressible  turbulent  boundary- layer  flow 
over  an  infinite  yawed  airfoil  has  been  examined  in  many  previous  in¬ 
vestigations;  e.  g. ,  see  Refs.  4,  7,  8,  9,  and  10.  The  important  finding 
from  these  studies  is  that  quite  small  departures  from  the  situation  of 
the  infinite  swept  wing  can  produce  appreciable  differences  in  boundary- 
layer  parameters.  This  remarkable  sensitivity  appears  to  be  directly 
traceable  to  small  spanwise  pressure  gradients  which  invalidate  the 
infinite  yawed  body  assumption  of  zero  spanwise  gradients.  The  pos¬ 
sibility  of  this  occurrence  must  be  kept  in  mind  when  one  compares 
theory  with  experiment  for  so-called  infinite  yawed  airfoil  flows. 

3.2.1  Comparison  with  Other  Analytical  Techniques 

Both  Nash  (Ref.  8)  and  Bradshaw  (Ref.  9)  have  developed  numerical 
calculation  techniques  for  the  three-dimensional  turbulent  boundary  lay¬ 
er  on  infinite  swept  wings  based  on  a  kinetic- energy- of- turbulence  ap¬ 
proach.  Reasonable  agreement  with  experiment  has  been  reported  for 
both  methods.  A  more  recent  investigation  by  Nash  and  Tseng  (Ref.  10) 
using  the  kinetic- energy-of- turbulence  approach  (in  conjunction  with 
the  assumption  that  the  turbulent  shear  stress  acts  in  the  direction  of 
the  mean  rate  of  strain)  is  concerned  with  high  Reynolds  number  incom¬ 
pressible  flow  over  a  35-deg  infinite  swept  wing  having  a  symmetrical, 

12 -per cent-thick  NACA  65iA012  airfoil  section  in  the  chordwise  direc¬ 
tion.  Calculations  were  performed  for  the  upper  surface  of  the  wing 
at  two  angles  of  incidence  -  corresponding  to  two-dimensional  lift  coef¬ 
ficients  of  0.  15  and  0.  50  -  and  for  two  values  of  the  free- stream  chord 
Reynolds  number,  namely  Re^c  =  6.  14  x  107,  which  corresponds  to  a 
full-scale  condition  for  a  large  transport  aircraft,  and  Re^c  =  0.  82  x 
107,  which  corresponds  to  a  typical  wind  tunnel  scale  for  an  aircraft 
model  or  to  full  scale  for  a  helicopter  rotor. 

In  order  to  assess  the  present  calculation  technique  using  an  eddy 
viscosity  approach  relative  to  the  kinetic- energy-of- turbulence  analysis 
by  Nash  and  Tseng  (Ref.  10),  turbulent  boundary- layer  calculations  have 
been  performed  for  the  upper  surface  of  the  above-described  NACA 
65jA012  airfoil  at  35-deg  sweep  with  a  sectional  lift  coefficient  of 
0.  50  and  a  free- stream  chord  Reynolds  number  of  Re0jc  =  6.14  x  107. 
The  inviscid  flow  field  is  determined  as  in  Section  2.  10  in  the  present 
report  (i.e.,  surface  pressure  distribution  is  derived  from  simple  sweep 
theory  using  the  two-dimensional,  potential-flow  velocity  distribution 
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tabulated  on  p.  370  of  Ref.  28).  A  plot  of  the  resulting  inviscid  surface 
velocity  distribution  may  be  found  in  Fig.  2  of  Nash  and  Tseng  (Ref.  10); 
this  figure  shows  an  adverse  pressure  gradient  over  almost  the  entire 
airfoil  with  only  the  leading- edge  region  in  a  favorable  pressure  gradi¬ 
ent.  In  order  to  be  consistent  with  the  analysis  by  Nash  and  Tseng  (Ref. 
10),  the  leading  edge  of  the  airfoil  is  herein  taken  to  be  turbulent. 

As  shown  in  Fig.  9,  agreement  between  the  present  infinite  yawed 
wing  calculation,  based  on  a  scalar  eddy  viscosity  model  of  turbulence, 
and  the  more  elaborate  kinetic- energy- of- turbulence  analysis  of  Nash 
and  Tseng  (Ref.  10)  is  excellent.  Recall  from  Section  2.  7  that  it  is^the 
chordwise  (x- direction)  component  of  the  displacement  thickness,  6X, 
which  has  physical  meaning  for  infinite  yawed  wing  flows,  which  ex¬ 
plains  the  choice  of  6*  for  comparison  in  Fig.  9.  Another  important 
point  is  that  the  present  analysis  technique  actually  integrates  the  gov¬ 
erning  turbulent  boundary- layer  equations  in  the  near- wall  region  using 
the  van  Driest  modification  of  the  eddy  viscosity;  Nash  and  Tseng  (Ref. 
10)  match  their  outer  numerical  solution  onto  an  inner-layer  solution 
which  takes  the  form  of  some  variation  of  the  universal  inner  law.  As 
discussed  by  Nash  in  Ref.  33,  a  generalized  version  of  Townsend's  in¬ 
ner  law  (Ref.  34)  is  used  with  modifications  which  permit  departures 
from  collateral  flow  through  the  inner  layer  as  well  as  small  deviations 
from  the  strict  linearity  of  the  variation  of  shear  stress  with  normal 
distance  from  the  wall.  It  is  believed  that  the  small  discrepancy  in 
both  chordwise  and  spanwise  skin-friction  coefficients  presented  in 
Fig.  9  may  be  a  function  of  the  assumed  inner  law  for  the  Nash  and 
Tseng  analysis. 

Further  comparison  of  the  present  scalar  eddy  viscosity  approach 
relative  to  the  Nash  and  Tseng  kinetic- energy-of-turbulence  analysis  is 
given  in  Fig.  10  relative  to  the  surface  streamline  direction,  us,  as  de¬ 
fined  in  Fig.  2.  Agreement  between  the  two  methods  is  excellent,  with 
the  choice  of  the  inner  law  as  discussed  in  the  previous  paragraph  being 
of  importance  relative  to  the  Nash  and  Tseng  results.  Note  from  Fig. 

10  that  the  three-dimensional  turbulent  boundary  layer  is  nearly  collat¬ 
eral  with  the  inviscid  outer- edge  flow  over  the  forward  half  of  the  wing; 
this  correspondence  is  caused  by  the  extreme  thinness  of  the  high ’Rey¬ 
nolds  number  turbulent  boundary  layer  in  this  region.  Over  the  rear  of 
the  wing,  as  the  boundary  layer  becomes  thicker  due  to  its  upstream 
growth,  the  adverse  pressure  gradient  causes  a  rapid  increase  in 
boundary- layer  crossflow,  which  is  reflected  as  an  increasing  differ¬ 
ence  between  the  surface  and  outer- edge  streamline  angles  in  Fig.  10. 
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The  most  favorable  comparisons  presented  above  indicate  that  the 
current  scalar  eddy  viscosity  model  of  three-dimensional  turbulent 
boundary- layer  flow  over  an  infinite  yawed  wing  yields  almost  identical 
results  to  the  three-dimensional  kinetic -energy- of- turbulence  analysis 
by  Nash  (Refs.  8  and  33)  as  applied  in  Nash  and  Tseng  (Ref.  10).  Hence 
it  must  be  concluded  that  the  present  scalar  eddy  viscosity  model  em¬ 
ploying  a  mixing-length  hypothesis  and  the  kinetic-energy-of -turbulence 
hypothesis  (Refs.  8,  9,  10,  and  33)  are  equally  satisfactory  for  numeri¬ 
cal  calculation  of  infinite  yawed  wing  turbulent  boundary-layer  flows 
under  adverse  pressure  gradient  situations  which  are  of  interest  for 
transport  aircraft  applications.  Similar  conclusions  have  been  reached 
by  Sivasegaram  and  Whitelaw  (Ref.  35),  who  compared  the  Patankar- 
Spalding  mixing-length  approach  (Ref.  17)  to  the  Bradshaw -Ferris  s- 
Atwell  kinetic-energy-of -turbulence  approach  (Ref.  36)  for  the  case 
of  two-dimensional  supersonic  turbulent  boundary  layers  in  favorable, 
zero,  and  adverse  pressure  gradient  flows  both  with  and  without  heat- 
transfer  effects . 

i 

3.2.2  Comparison  with  Experiment 

Probably  the  most  complete  and  best- documented  set  of  experi¬ 
mental  data  for  subsonic  flow  over  an  infinite  yawed  airfoil  is  that  re¬ 
ported  by  Altman  and  Hayter  (Ref.  37)  some  twenty  years  ago.  This 
experimental  investigation  was  conducted  in  the  Ames  Aeronautical 
Laboratory  7-  by  10-ft  subsonic  wind  tunnel  using  both  swept  and  un¬ 
swept  NACA  63j-012  section  airfoils.  The  unswept  wing,  having  a 
chord  of  4  ft,  was  mounted  vertically  and  spanned  the  7- ft  dimension 
of  the  wind  tunnel  (see  Fig.  2  of  Ref.  37  for  clarification).  The  swept 
wing,  having  a  sweep  angle  of  45  deg  and  a  chord  of  2.  5  ft  perpendicular 
to  the  leading  edge,  was  mounted  horizontally  and  spanned  the  10-ft 
dimension  of  the  wind  tunnel  (see  Fig.  3  of  Ref.  37  for  clarification). 
Detailed  boundary- layer  velocity  profiles  were  measured  along  the  mid¬ 
span  section  at  30,  40,  50,  60,  70,  80,  and  90  percent  of  the  stream- 
wise  chord  on  both  wings.  In  addition,  the  flow  direction  within  the 
boundary  layer  was  determined  for  the  swept  wing.  These  measure¬ 
ments  were  made  for  section  lift  coefficients  of  0,  0.  32,  0.  46,  0.  74, 
and  1.  00  based  on  the  component  of  velocity  normal  to  the  leading  edge. 
In  conjunction  with  the  above- described  profile  measurements,  the  dis¬ 
tribution  of  surface  pressure  over  both  wings  was  also  measured;  for 
the  swept  wing  the  spanwise  distribution  of  surface  pressure  was  meas¬ 
ured  along  constant- chord  lines  at  5,  15,  30,  50,  and  80  percent  of  the 
chord.  The  tests  of  the  swept  wing  were  conducted  at  a  dynamic  pres¬ 
sure  of  150  psf  and  those  on  the  unswept  wing  at  a  dynamic  pressure 
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of  29.  5  psf.  The  free- stream  Reynolds  number  based  on  the  airfoil 
chord  and  the  component  of  free- stream  velocity  normal  to  the  leading 
edge  was  approximately  3.  8  x  10®  for  both  wings. 

A  comparison  of  experimental  and  theoretical  chordwise  pressure 
distributions  on  the  upper  surface  of  the  NACA  63j-012  airfoil  at  45- 
deg  sweep  in  presented  in  Fig.  11  for  two  values  of  the  sectional  lift 
coefficient,  Clu,  namely  0.  0  and  0.  74.  The  theoretical  distribution  is 

based  upon  application  of  classical  sweep  theory  (as  previously  discussed 
in  Section  2.  10)  using  the  surface  velocity  distribution  tabulated  in  Table 
I,  which  is  taken  directly  from  p.  337  of  Abbott  and  von  Doenhoff  (Ref. 
28).  As  is  obvious  from  Fig.  11,  classical  sweep  theory  is  in  excellent 
agreement  with  experiment  over  the  entire  airfoil  for  both  values  of  the 
sectional  lift  coefficient.  Further  comparisons  of  sweep  theory  relative 
to  experiment  for  the  airfoil  of  current  interest  may  be  found  in  Ref.  38. 

Before  examining  the  three-dimensional  turbulent  boundary- layer 
results  for  the  swept  airfoil,  it  would  be  well  to  first  ascertain  the  va¬ 
lidity  of  the  present  analysis  under  two-dimensional  conditions  of  zero 
sweep  and  zero  lift.  Presented  in  Fig.  12  is  a  comparison  of  calculated 
chordwise  turbulent  boundary- layer  parameters  (displacement  thick¬ 
ness  6*,  momentum  thickness  0mx,  ancl  shape  factor  Hfx)  relative  to 
experimental  measurements  for  zero  sweep  and  zero  yaw  as  given  in 
Fig.  7  of  Ref.  37.  Instantaneous  transition  from  laminar  to  turbulent 
flow  is  taken  to  occur  at  20  percent  of  the  chord,  which  corresponds  to 
the  location  of  artifically  induced  (roughness)  transition  in  the  experi¬ 
ments  of  Ref.  37  under  zero-lift  conditions.  Agreement  between  the 
present  analysis  technique  and  experiment  is  excellent  for  this  two- 
dimensional  condition.  Further  illustration  of  the  applicability  and 
validity  of  the  present  turbulent  boundary-layer  model  may  be  seen  in 
Fig.  13,  which  compares  a  calculated  turbulent  boundary- layer  velo¬ 
city  profile  with  experimental  probe  measurements  (see  Fig.  6  of  Ref. 

37)  for  the  50-percent  chord  location  on  the  zero  sweep  -  zero  yaw 
airfoil. 

Turn  now  to  three-dimensional  flow  situations.  Figure  14  shows 
the  calculated  chordwise  (x- direction)  turbulent  boundary- layer  param¬ 
eters  (displacement  thickness  6*,  momentum  thickness  0m,x»  and 
shape  factor  Hf  x)  relative  to  experimental  measurements  as  given  in 
Fig.  7  of  Ref.  37  for  the  case  of  45- deg  sweep  and  zero  lift.  Instantan¬ 
eous  transition  from  laminar  to  turbulent  flow  is  taken  to  occur  at  20- 
percent  chord,  which  corresponds  to  the  location  of  artifically  induced 
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(roughness)  transition  in  the  experiments  of  Ref.  37  under  zero-lift 
conditions.  Agreement  between  the  present  analysis  technique  and  ex¬ 
periment  is  excellent  for  this  yawed  wing  at  zero-lift  flow.  A  com¬ 
parison  of  calculated  chordwise  (u/Ue)  and  spanwise  (w/We)  turbulent 
boundary- layer  velocity  profiles  at  the  50-percent  chord  location  is 
given  in  Fig.  15  relative  to  experimental  measurements  reported  in 
Fig.  6  of  Ref.  37  for  the  45- deg  sweep  condition  at  zero  lift.  Note 
from  Fig.  15  that  the  chordwise  and  spanwise  velocity  profiles  are  al¬ 
most  identical,  which  is  a  direct  consequence  of  the  almost  constant 
surface  pressure  distribution  between  about  10-percent  chord  and  the 
50-percent  chord  location  of  present  interest  (see  the  Cj_,  =  0.  0  curve 
given  in  Fig.  11  for  clarification), 

A  two-dimensional  turbulent  calculation  in  the  chordwise  direction 
would  predict  the  chordwise  components  of  the  boundary-layer  flow  cor¬ 
rectly  if  the  so-called  principle  of  independence  held,  as  it  does  in  the 
case  of  incompressible  laminar  boundary- layer  flow  over  infinite  yawed 
cylinders  and  wings.  According  to  the  principle  of  independence,  the 
shear- stress  component  in  a  particular  direction  depends  only  on  the 
resolved  velocity  profile  in  that  direction.  The  failure  of  the  principle 
of  independence  in' the  turbulent  case  has  received  wide  coverage  in  the 
literature,  and  it  is  now  well  known  that  a  chordwise  two-dimensional 
calculation  does  not  predict  the  correct  behavior  of  the  turbulent  bound¬ 
ary  layer  on  a  yawed  flat  plate.  (See  pp.  54  through  56  of  Ref.  39  for 
a  good  discussion  of  this  point.  )  As  can  be  seen  from  a  review  of  Sec¬ 
tions  2.  3,  2.  4,  and  2.  5  in  the  present  report,  the  scalar  eddy  viscosity 
model  of  turbulence  as  used  herein  does  not  admit  the  principle  of  in¬ 
dependence  since,  for  example,  the  chordwise  shear  stress,  rx,  as 
given  by  Eq.  (15)  is  a  function  of  both  the  chordwise  and  spanwise  ve¬ 
locity  profiles  through  the  eddy  viscosity  relationship  given  by  Eqs. 

(17)  and  (18);  the  same  is  true  for  the  spanwise  shear  stress,  rz, 
given  by  Eq.  (16).  Practical  examples  illustrating  the  effects  of  sweep 
on  the  present  airfoil  at  zero-lift  flow  situation  are  given  in  Figs.  16 
and  17.  Note  from  Fig.  16  that  the  influence  of  sweep  (and  hence  three- 
dimensionality)  is  to  increase  the  magnitude  of  the  chordwise  turbulent 
boundary-layer  properties  at  a  given  chordwise  location.  Further  ob¬ 
serve  from  Fig.  17  that  the  influence  of  sweep  is  to  increase  the  chord- 
wise  turbulent  boundary- layer  thickness  at  the  50-percent  chord  loca¬ 
tion.  Similar  findings  have  been  reported  by  Ashkenas  and  Riddell 
(Ref.  40)  for  turbulent  boundary- layer  flow  over  yawed  and  unyawed 
flat  plates. 


30 


AEDC-TR-73-112 


As  discussed  previously  in  this  section,  the  laminar  boundary 
layer  was  tripped  to  turbulent  flow  at  20-percent  chord  location  for 
both  the  yawed  and  the  zero-yaw  flows.  This  means  that  natural  tran¬ 
sition  did  not  occur  under  the  experimental  flow  conditions  of  present 
interest.  Recall  from  Section  2.  8  that  natural  boundary- layer  transi¬ 
tion  on  yawed  wings  appears  to  be  governed  by  a  three-dimensional 
crossflow  instability  phenomenon  which  can  be  correlated  with  a  so- 
called  maximum  local  crossflow  Reynolds  number,  X,  defined  by  Eq. 
(97);  all  of  the  work  to  date  on  crossflow- dominated  natural  transition 
indicates  that 


X  <  100 — ►Laminar  Boundary  Layer 

100  ^  X  ^  200 — ►Vortex  Formation  and  Transitional 

Boundary  Layer 

X  >  200 — ►Turbulent  Boundary  Layer 

As  shown  in  Fig.  18  the  calculated  crossflow  Reynolds  number,  X,  dis¬ 
tribution  for  the  present  flow  situation  (NACA  63j-012  airfoil  at  45-deg 
sweep  with  Re^c  =  5.  37  x  10^)  attains  a  maximum  value  of  approxi¬ 
mately  120  at  a  chordwise  location  of  X/C  =  0.  27.  As  shown  in  the  a- 
bove  discussion,  crossflow- dominated  natural  transition  should  not  oc¬ 
cur  on  the  present  wing;  this  is  in  agreement  with  the  experimental  ob¬ 
servations  of  Ref.  37.  It  should  be  further  noted  from  Fig.  18  that  the 
physical  chordwise  location  of  the  artificial  roughness  used  to  promote 
boundary- layer  transition  coincides  fairly  closely  with  the  chordwise 
location  of  maximum  crossflow  Reynolds  number,  X,  where  it  is  to  be 
expected  that  the  three-dimensional  laminar  boundary  layer  is  most 
unstable  with  respect  to  external  disturbances.  The  author  believes 
that,  generally,  boundary- layer  trips  on  yawed  wings  should  be  placed 
to  correspond  with  the  chordwise  location  of  maximum  crossflow  Rey¬ 
nolds  number,  X,  in  order  to  promote  effective  tripping. 

All  of  the  above  results  and  discussion  have  been  for  the  zero-lift 
airfoil  under  both  yawed  and  unyawed  conditions.  For  transport  air¬ 
craft  applications  the  important  area  of  interest  to  the  aerodynamic ist 
is  the  upper  surface  of  swept  lifting  airfoils  relative  to  adverse  pres¬ 
sure  gradient  effects  on  the  three-dimensional  turbulent  boundary  layer. 
The  present  45-deg  sweep  NACA  63^-012  airfoil  at  angle  of  attack  such 
that  the  sectional  lift  coefficient,  Cj^,  equals  0.  74  is  an  excellent  ex¬ 
ample  of  this  type  of  flow.  Reference  to  Fig.  11  presented  previously 
shows  that  the  pressure  gradient  over  the  entire  upper  surface  of  the 
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airfoil  is  adverse  with  a  particularly  large  gradient  in  the  leading- edge 
region.  Comparisons  of  calculated  and  experimental  turbulent  boundary- 
layer  parameters  (displacement  thickness,  6*,  momentum  thickness, 

0m,  and  shape  factor,  Hf)  for  both  the  chordwise  (x- direction)  and 
spanwise  (z-direction)  components  of  flow  over  the  upper  surface  of 
the  airfoil  are  given  in  Fig.  19.  Boundary-layer  transition  fpom  lami¬ 
nar  to  turbulent  flow  is  taken  to  occur  near  the  leading  edge  of  the  air¬ 
foil  in  accord  with  the  experimental  results  of  Ref.  37.  As  can  be  seen 
from  Fig.  19,  agreement  between  the  present  analysis  technique  and  ex¬ 
periment  is,  in  general,  good  over  the  entire  airfoil  with  the  largest 
discrepancy  in  the  chordwise  displacement  thickness  near  the  trailing 
edge.  The  chordwise  and  spanwise  turbulent  boundary- layer  velocity 
profiles  at  the  50-percent  chord  location  are  presented  in  Fig.  20  for 
the  current  yawed  lifting  airfoil  condition.  The  calculated  spanwise 
velocity  profile,  w/We,  is  in  excellent  agreement  with  experiment; 
the  calculated  chordwise  velocity  profile,  u/Ue,  is  somewhat  different 
from  the  experimental  profile  in  the  near-wall  region  but  is  in  excellent 
agreement  near  the  outer  edge  of  the  boundary  layer.  Adverse  pres¬ 
sure  gradient  effects  on  both  the  chordwise  and  spanwise  turbulent 
boundary- layer  velocity  profiles  become  apparent  upon  examination  of 
Figs.  21  and  22,  which  are  for  the  50-percent  chord  and  trailing- edge 
locations,  respectively,  on  the  upper  surface  of  the  present  lifting 
swept  airfoil.  The  spanwise  profile  remains  typical  of  a  zero  pressure 
gradient  turbulent  flow,  whereas  the  chordwise  profile  clearly  illus¬ 
trates  the  retarding  effects  of  the  chordwise  adverse  pressure  gradient. 
Corresponding  calculated  eddy  viscosity  distributions  across  the  turbu¬ 
lent  boundary  layer  at  the  above  two  locations  (50-  and  100-percent  chord) 
are  given  in  Fig.  23.  Also  indicated  on  this  figure  are  inner  and  outer 
regions  which  correspond  to  the  particular  form  of  equation  used  for  the 
eddy  viscosity  relationship,  namely  Eq.  (28)  for  the  inner  region  calcu¬ 
lation  and  Eq.  (29)  for  the  outer  region  calculation.  As  is  obvious  from 
Fig.  23,  the  turbulence  levels  in  the  outer  region  of  the  boundary  layer 
are  quite  large. 

Comparisons  of  calculated  streamline  direction,  w,  within  the 
three-dimensional  turbulent  boundary  layer  at  various  chordwise  loca¬ 
tions  on  the  upper  surface  of  the  above- discussed  swept  lifting  airfoil 
are  given  in  Fig.  24  relative  to  experimental  measurements  (via  a  di¬ 
rectional  rake)  reported  in  Fig.  16  of  Ref.  37.  Agreement  between  the 
present  analysis  technique  and  experiment  is  reasonably  good,  with  the 
largest  discrepancy  in  the  inner  region  of  the  boundary  layer  near  the 
wall.  This  behavior  can  be  traced  directly  to  the  difference  between 
calculated  and  experimental  chordwise  (x-direction)  velocity 
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distributions  near  the  wall  as  discussed  previously  in  connection  with 
Fig.  20.  The  physical  reason  for  this  near-wall  disagreement  may  well 
be  the  small  spanwise  pressure  gradient  over  the  swept  lifting  airfoil 
under  experimental  conditions  as  discussed  in  Ref.  37.  Recall  from 
the  beginning  of  Section  3.  2  that  small  spanwise  pressure  gradients  can 
cause  large  departures  from  the  true  infinite  swept-wing  situation,  re¬ 
sulting  in  appreciable  differences  relative  to  boundary-layer  parameters. 


3.3  INFINITE  EXTENT  YAWED  AIRFOIL  UNDER  TRANSONIC  CONDITIONS 

To  the  author's  knowledge  the  case  of  transonic  compressible 
boundary- layer  flow  over  an  infinite  yawed  airfoil  has  not  been  consid¬ 
ered  in  any  previous  investigation,  either  analytical  or  experimental. 
The  recent  report  by  Nash  and  Scruggs  (Ref.  3)  describes  the  calcula¬ 
tion  of  the  three-dimensional  compressible  turbulent  boundary  layer  on 
the  supercritical  finite  swept  wing  of  the  NASA  modified  F8  transonic 
research  aircraft  using  the  kinetic- energy- of- turbulence  method  of 
Nash  (Refs.  2,  8,  and  33)  extended  to  compressible  flow  via  the  Crocco 
relationship  for  temperature  (which  restricts  the  analysis  to  adiabatic- 
wall  flows).  No  comparisons  with  experiment,  either  wind  tunnel  or 
flight,  are  presented  in  Ref.  3  since  such  data  are  classified.  This 
particular  area  of  transonic  research  is  currently  undergoing  much  ef¬ 
fort  for  transonic  transport  applications  and  hence  will  not  be  further 
considered  in  the  present  report.  Instead,  the  following  treatment  of 
hot-wall  transonic  boundary-layer  flow  will  be  pursued  because  of  its 
direct  applicability  to  the  design  of  ground  test  experiments  for  cor¬ 
rect  simulation  of  aerodynamic  phenomena  to  be  encountered  during 
transonic  and  subsonic  flight  of  lifting  entry  vehicles  such  as  the  NASA 
Space  Shuttle. 

With  the  advent  of  maneuverable  reentry  vehicles  such  as  the  Space 
Shuttle,  a  new  problem  involving  heat  transfer  from  hot  aerodynamic 
surfaces  has  evolved.  During  entry  into  the  atmosphere,  the  Shuttle 
will  be  subjected  to  large  heat  loads  which  depend  upon  the  geometric 
configuration  as  well  as  the  attitude  of  the  vehicle  during  reentry. 

Some  of  this  entry  heat  load  will  inevitable  soak  into  the  structure  and 
aerodynamic  surfaces.  After  the  hypersonic  reentry  phase  of  the 
flight,  the  Shuttle  will  slow  to  subsonic  speeds  and  maneuver  for  a  con¬ 
ventional  aircraft  landing.  It  has  been  estimated  (see  Ref.  41)  that  the 
Shuttle  wings  may  reach  soak  temperatures  on  the  order  of  1000°R, 
whidh  is  some  2.  5  times  the  local  free-stream  stagnation  temperature 
under  subsonic  flight  conditions.  Thus,  during  the  transonic 
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maneuvering  and  subsonic  landing  phase  of  Shuttle  reentry,  the  heated 
aerodynamic  surfaces  will  be  (literally)  injecting  heat  into  a  high  Rey¬ 
nolds  number  turbulent  boundary  layer.  The  important  question  rela¬ 
tive  to  transonic  testing  of  Shuttle  vehicles  becomes,  "What  are  the 
aerodynamic  effects  of  high  rates  of  heat  transfer  from  the  surface  to 
a  transonic,  turbulent  boundary  layer,  and  must  we  simulate  this  hot- 
wall  environment  in  ground  testing  of  Shuttle  configurations?"  Of  par¬ 
ticular  interest  is  the  effect  of  a  hot  wall  on  the  three-dimensional 
transonic  turbulent  boundary- layer  which  will  be  herein  considered 
through  examination  of  infinite  yawed  airfoil  flows.  To  the  author's 
knowledge,  the  present  analysis  technique  is  the  only  available  method 
which  allows  calculation  of  infinite  yawed  airfoil  flows  including  flow 
compressibility  and  surface  heat-transfer  effects.  In  this  connection 
the  current  method,  slightly  modified  with  respect  to  the  inviscid  flow- 
field  input  to  the  boundary  layer,  has  been  applied  under  supersonic 
and  hypersonic  conditions  with  excellent  results.  (See  Ref.  25  for  com¬ 
plete  documentation  relative  to  Space  Shuttle  configurations.) 

In  order  to  provide  some  insight  into  the  above-defined  problem 
area,  to  which  little  or  no  attention  has  been  given  in  previous  studies, 
infinite  yawed  wing  calculations  based  on  the  present  analysis  techni¬ 
que  have  been  performed  for  the  upper  surface  of  the  Douglas  Aircraft 
Company  (DAC)  30- deg- sweep  airfoil  given  in  Fig.  15  of  Ref.  1.  Free- 
stream  conditions  used  in  the  calculations  are  as  follows: 

M.  =  0.  815 

Re<BjC  =  1.  05  x  108 

T0ja>  =  489°R,  T.  =  432°R 

po,®  =  325  psia,  P®  =  210  Psia 

C  =  1.  0  ft 

These  values  are  representative  of  the  proposed  AEDC  high  Reynolds 
number  transonic  tunnel  (HIRT)  as  described  in  Ref.  42..  The  above- 
defined  free- stream  Mach  number/ chord  Reynolds  number  combination 
is  very  close  to  actual  flight  values  quoted  in  Fig.  3  of  Ref.  43  for  the 
McDonnell  Douglas  Astronautics  Company  (MDAC)  Shuttle  Orbiter  at 
an  altitude  of  42,  500  ft.  In  all  of  the  calculated  results  to  be  presented 
below,  hot-wall  (surface  temperature  of  1000°R)  results  typical  of  ac¬ 
tual  Shuttle  entry  conditions  will  be  compared  with  adiabatic-wall  results 
since  an  adiabaticnvall  condition 'represents  conventional  transonic  test¬ 
ing  in  continuous  wind  tunnels  such  as  the  AEDC  Propulsion  Wind  Tun¬ 
nel  (16T)  and  Aerodynamic  Wind  Tunnel  (4T). 
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Presented  in  Fig.  25  are  the  inviscid  flow  parameters  over  the  up¬ 
per  surface  of  the  DAC  airfoil.  As  can  be  seen  from  the  figure,  the 
chordwise  (x-direction)  inviscid  flow  is  supercritical  (local  x-direction 
Mach  number  exceeds  the  sonic  of  critical  value),  whereas  the  span- 
wise  (z-direction)  inviscid  flow  remains  subsonic  throughout.  For  the 
present  study  the  chordwise  inviscid  supercritical  flow  is  not  allowed 
to  terminate  via  a  shock  wave  but  is  forced  to  recompress  smoothly 
over  the  aft  section  of  the  airfoil  with  a  continually  increasing  adverse 
pressure  gradient;  this  is  done  in  order  to  assess  the  effects  of  con¬ 
tinuous  adverse  pressure  gradient  on  three-dimensional  transonic  tur¬ 
bulent  boundary- layer  flow.  Note  in  Fig.  25  that  the  airfoil  chord 
length  is  one  foot,  which  is  representative  of  allowable  model  size  in 
the  proposed  AEDC  H1RT  facility. 

Figure  26  presents  the  calculated  distribution  for  the  ratio  of  the 
adiabatic -wall  temperature  to  the  free- stream  total  temperature.  Note 
that  this  ratio  ranges  from  0.  975  to  0.  995  with  a  distribution  which  fol¬ 
lows  in  trend  the  surface  pressure  distribution  of  the  previous  figure. 

The  hot-wall  condition  (Tw  =  1000°R)  yields  a  wall  temperature  ratio 
(Tw/T0  a3)  of  2.  043,  which  remains  constant  over  the  entire  airfoil. 

Figures  27  through  31  illustrate  the  results  of  the  hot-  (1000°R) 
versus  adiabatic-wall  calculations  of  present  interest.  As  can  be  seen 
from  Fig.  27,  three-dimensional  turbulent  boundary-layer  separation 
(defined  as  the  chordwise  surface  location  where  the  surface  flow  angle, 
<os,  is  identically  zero  following  Refs.  39  and  44)  occurs  at  93-percent 
chord  for  the  hot-wall  case  and  99-percent  chord  for  the  adiabatic -wall 
condition;  note  that  over  the  entire  aft  half  of  the  airfoil  the  adiabatic- 
wall  surface  flow  angle  is  always  larger  than  the  corresponding  x- 
location  hot- wall  value.  The  chordwise  (x-direction)  and  spanwise  (z- 
direction)  components  of  the  three-dimensional  turbulent  skin-friction 
coefficient  are  presented  in  Fig.  28,  with  the  important  point  being 
that  the  hot  wall  reduces  the  local  x-direction  skin  friction  by  some  20 
to  40  percent  over  the  corresponding  location  x-direction  adiabatic  wall 
value.  Similarly,  the  z-direction  hot- wall  skin  friction  is  reduced 
some  20  to  25  percent  over  the  entire  airfoil  as  compared  with  the  cor¬ 
responding  location  z-direction  adiabatic-wall  value.  Needless  to  say, 
this  significant  decrease  in  local  skin  friction  due  to  the  heated  wall 
can  result  in  appreciable  reduction  of  integrated  skin-friction  drag  on 
the  airfoil.  The  corresponding  calculated  three-dimensional  turbulent 
boundary- layer  parameters  (Hf,  6*,  0m)  for  both  the  x-  and  z-directions 
are  presented  in  Figs.  29  through  31.  Note  from  Fig.  29  the  sizable  in¬ 
crease  in  the  shape  factors  Hf<x  and  Hf  z  caused  by  the  heated  wall  as 
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compared  with  the  adiabatic  wall.  The  reason  for  this  behavior  can  be 
seen  from  Figs.  30  and  31  where,  in  general,  the  hot  wall  increases 
the  displacement  thicknesses  6*  and  6Z  over  the  entire  airfoil  but  de¬ 
creases  the  momentum  thicknesses  0mjX  and  ®rn  z  over  the  majority 
of  the  airfoil.  Recalling  from  Section  2.  7  that  it  is  the  x- direction  dis¬ 
placement  thickness,  6^,  which  has  physical  meaning  relative  to  infi¬ 
nite  yawed  wing  boundary-layer  interaction  with  the  inviscid  flow.  Fig. 

30  shows  that  the  hot  wall  increases  the  magnitude  of  6*  some  50  per¬ 
cent  over  the  entire  airfoil  as  compared  with  the  corresponding  x- 
location  adiabatic-wall  value. 

The  calculated  surface  heat-transfer  rate  as  reflected  through  the 
Stanton  number,  St,,,,  defined  by  Eq.  (89)  is  given  in  Fig.  32  for  the 
present  hot-wall  condition.  As  shown  on  the  figure,  the  local  heat  flux, 
qw,  is  positive  (denoting  heat  flow  from  the  airfoil  to  the  flow)  while 
the  wall  temperature  is  greater  than  the  stagnation  temperature  (Tw/ 

T0  m  =  2.043)  which  together  combine  to  yield  a  positive  value  for  the  Stan¬ 
ton  number.  Note  the  high  rate  of  heat  transfer  in  the  leading- edge 
region  with  essentially  a  constant  heat-transfer  rate  over  the  mid¬ 
section  of  the  airfoil  and  a  decrease  in  the  heating  rate  as  the  trailing 
edge  is  approached.  Recall  that  three-dimensional  boundary-layer 
separation  occurred  at  93-percent  chord. 

Hot- wall  effects  on  the  three-dimensional  turbulent  boundary-layer 
velocity  and  temperature  profiles  are  illustrated  in  Figs.  33  and  34  rel¬ 
ative  to  the  adiabatic- wall  results  at  a  common  chordwise  location, 
namely  90- percent  chord.  This  particular  location  is  chosen  to  permit 
examination  of  turbulent  boundary- layer  structure  having  an  upstream 
adverse  pressure  gradient  history.  As  is  obvious  from  Fig.  33,  the 
hot  wall  tends  to  retard  the  chordwise  velocity  profile,  u/Ue,  more 
than  does  the  adiabatic  wall;  this  results  in  earlier  boundary- layer 
separation  under  hot- wall  conditions,  as  noted  earlier.  The  spanwise 
velocity  profile,  w/We,  is  not  strongly  affected  by  the  hot  wall  except 
in  the  near-wall  region,  where  the  van  Driest  damping  term  given  by 
Eq.  (27)  becomes  important.  The  most  important  hot-wall  effects  on 
turbulent  boundary- layer  structure  are  seen  in  Fig.  34  with  respect  to 
both  the  static  temperature  ratio,  T/Te,  and  stagnation  temperature 
ratio,  T0/T0j00.  Due  to  the  constancy  of  static  pressure  across  the 
boundary  layer,  the  static  temperature  ratio,  T/Te,  is  the  inverse  of 
the  static  density  ratio,  p/pe,  through  the  perfect  gas  equation  of  state, 
p  =  pRT.  Hence  the  hot-wall  case  under  present  examination  results 
in  appreciable  density  variation  across  the  turbulent  boundary  layer  as 
compared  with  the  adiabatic-wall  condition.  This  density  variation 
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becomes  directly  reflected  in  the  eddy  viscosity  distribution  across  the 
boundary  layer  as  can  be  seen  from  Eqs.  (28)  and  (29). 

All  of  the  above-presented  results  clearly  indicate  that  hot-wall 
effects  can  indeed  have  a  significant  influence  on  transonic  aerodynam¬ 
ics.  The  recent  experimental  investigation  of  Ref.  41  conducted  in  the 
Texas  A  &  M  University  7-  by  10-ft  subsonic  wind  tunnel  using  an  elec¬ 
trically  heated  NACA  0012-64  airfoil  reveals  that  both  the  maximum 
lift  coefficient,  CjLmax'  anc*  t*ie  sta^  angle  of  the  airfoil,  astaip  are 

reduced  as  the  wing-to-free -stream  temperature  ratio,  TW/T0D,  is  in¬ 
creased  under  subsonic  conditions;  e.g.,: 


Tw/T„ 

CI 

■‘-'max 

“stall 

1.0 

1.  6 

17.  0 

2.  0 

1.  2 

— 

2.2 

— 

12.  5 

Hence,  simulation  of  the  hot  wall/cold  free- stream  environment  using 
electrically  heated  models  appears  necessary  in  ground  testing  under 
subsonic  and  transonic  conditions  for  Space  Shuttle  applications,  espe¬ 
cially  for  aerodynamic  drag  and  stall  characteristics.  More  work,  both 
analytical  and  experimental,  is  definitely  needed  to  fully  explore  and  un¬ 
derstand  this  relatively  new  area  of  hot- wall  subsonic  and  transonic 
aerodynamics  for  application  to  lifting  bodies  such  as  the  Space  Shuttle. 
Especially  desirable  for  future  studies  are  turbulent  boundary- layer 
structure  measurements  in  two-  and  three-dimensional  subsonic  and 
transonic  flows  under  controlled  hot- wall  conditions. 


SECTION  IV 

CONCLUDING  SUMMARY 


The  present  report  has  documented  a  three-dimensional  compress¬ 
ible  turbulent  boundary- layer  analysis  applicable  to  yawed  airfoil 
flows  of  infinite  extent  under  subsonic  and  transonic  conditions.  The 
boundary- layer  analysis  has  been  based  on  implicit  finite-difference 
integration  of  the  governing  infinite  yawed  body  boundary-layer  equa¬ 
tions  for  turbulent  flow  using  a  scalar  eddy  viscosity  model  of  three- 
dimensional  turbulence.  Comparisons  with  other  analysis  techniques 
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as  well  as  experimental  .measurements  under  subsonic  wind  tunnel  con¬ 
ditions  have  been  presented  to  establish  and  ascertain  the  basic  validity 
and  applicability  of  the  current  technique  for  infinite  yawed- wing  flows. 

Ground  testing  of  Space  Shuttle  configurations  under  continuous 
transonic  flow  conditions  with  an  adiabatic  wall  may  not  be  totally  ap¬ 
plicable  to  actual  Shuttle  entry  where  the  wing  surface  temperature 
may  reach  soak  values  on  the  order  of  twice  the  free- stream  stagna¬ 
tion  temperature  because  of  the  hypersonic  high  heating  phase  of  the 
reentry  trajectory.  Numerical  calculations  presented  in  this  report 
indicate  an  appreciable  reduction  of  integrated  skin-friction  drag  on  a 
hot  wall  (Tw/T0  „  «  2)  airfoil  relative  to  the  same  airfoil  under  adia¬ 
batic  wall  conditions.  Airfoil  static  stall  angle  and  maximum  lift  coef¬ 
ficient.  as  well  as  turbulent  boundary-layer  parameters,  are  strongly 
affected  by  hot-wall  relative  to  adiabatic- wall  conditions.  Hence,  sim¬ 
ulation  of  the  hot  wall/cold  free- stream  environment  using  electrically 
heated  models  appears  necessary  in  ground  testing  under  continuous 
subsonic  and  transonic  wind  tunnel  conditions  for  Space  Shuttle  appli¬ 
cations,  especially  for  aerodynamic  drag  and  stall  characteristics 
determination. 

The  present  work  makes  no  attempt  to  address  the  problem  of 
viscous- inviscid  interactions  at  subsonic  or  transonic  speeds.  (The 
term  viscous- inviscid  interaction  is  taken  here  to  embrace  all  flow 
situations  in  which  the  development  of  the  boundary  layer  and  wake  has 
a  significant  effect  on  the  pressure  field. )  An  excellent  discussion  of 
this  topic  for  the  case  of  adiabatic -wall  transonic  flows  may  be  found 
in  Ref.  45. 
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Infinite-Extent  Yawed  Body 

W  -  W  ■ 
e  oo  oz 


Stagnation  or  Attachment  Line  (x  ■  0) 


(Free-Stream) 


x,  u;  Chordwise  Direction 


z,  w;  Spanwise  Direction 


Fig.  1  Infinite  Extent  Yawed  Body  Geometry  and  Nomenclature 


Limiting  Surface' 
Streamline  — ■ 


Chordwise  Flow  Separation  Line 
(where  uc  s  0) 


In  viscid  Edge 
Streamline 


Trailing  Edge 


a.  Plan  View  Illustrating  Three-Dimensional  Boundary-Layer 
Separation  on  an  Infinite  Yawed  Wing 


h.  Chordwise  Section  of  Wing 
Fig.  2  Infinite  Yawed  Wing  Geometry  and  Nomenclature 
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Fig.  3  Schematic  of  Stagnation  or  Attachment  Line  Flow  over 
an  Infinite  Extent  Yawed  Circular  Cylinder 
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Fig.  4  Coordinate  System  for  Flow  over  an  Infinite  Yawed  Circular  Cylinder 
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Twisted  Profile 


Plane  Tangential 
to  Outer  Flow 
Streamline 

Tangential 
Component ,  ug  ^ 


Plane  Normal  to  the  Outer 
Flow  Streamline ,  i.e.,  the 
Crossflow  Plane 


Crossflow  Component,  w  . 

s  t 

Profile  Inflection  Point 


si,  max 


Body 

Surface 


'si 


Fig.  5  Three-Dimensional  Boundary-Layer  Velocity  Profiles  in  Streamline  Coordinates 
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Stagnation  or  Attachment-Line  Flow 
on  an  Infinite  Yawed  Airfoil 

■  Present  Boundary-Layer  Calculations 

O  ▼  Experimental  Data  from  Fig.  4 


Fig.  6  Variation  of  Stagnation  or  Attachment  Line  Rea  e  m  f  with  C* 
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Fig.  8  Variation  of  Stagnation  or  Attachment  Line  C*  with  C# 
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foil  at  35-deg  Sweep 
07,  Adiabatic  Wall, 

t  Boundary-Layer  Theory 

>f-Turbulence  Analysis 
g  (Ref.  10) 


NACA  63^-012  Airfoil  at  45. 0-deg  Sweep 

Classical  Sweep  Theory  Using  Surface  Velocity 
Distribution  from  Table  i 


o  □  Experimental  Data  from  Figs.  5  and  10  in  Ref.  37 


j _ i _ I _ i _ i _ I _ i _ L 

10  20  30  40  50  60  70  80 

Percent  X/C 


Fig.  11  Surface  Pressure  Distribution  on  an  Infinite  Yawed  Airfoil  under  Subsonic  Conditions 
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Present  Turbulent  Boundary-Layer 
Theory 

O  Experimental  Data  from  Fig.  7  in 
0.004  r-  Ref.  37 


Percent  X/C 

Fig.  12  Turbulent  Boundary-Layer  Parameters  on  a  Two-Dimensional 
Airfoil  at  Zero  Lift 
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Present  Turbulent  Boundary -Layer  Theory 

Experimental  Data  from  Fig.  9  in 
Ref.  37 


Percent  X/C 

Fig.  14  Turbulent  Boundary- Layer  Parameters  on  an  Infinite  Yawed 
Airfoil  at  Zero  Lift 
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Fig.  15  Turbulent  Boundary-Layer  Velocity  Profiles  on  an  Infinite 
Yawed  Airfoil  at  Zero  Lift 
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0.005 


0.004 


0.003 

c r 

0.002 


0.001 


NACA  631-012  Airfoil 

■  Present  Turbulent  Boundary-Layer 
Theory  at  45.0-deg  Sweep  from  Fig.  14 

•———Present  Turbulent  Boundary-Layer 

Theory  at  0.0-deg  Sweep  from  Fig.  12 


1.6 


30  40  50  60  70  80  90  100 


Percent  X/C 

Fig.  16  Comparison  of  Two-Dimensional  and  Infinite  Yawed  Airfoil 
Boundary-Layer  Parameters 
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Fig.  18  Crossflow  Reynolds  Number  Distribution  on  an  Infinite  Yawed  Airfoil  at  Zero  Lift 
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NACA  63^-012  Airfoil  at  45. 0-deg  Sweep 
Reoo,  C  5.4x1^  Adiabatic  Wall,  Cl  ■  0. 74 

— —  Present  Turbulent  Boundary-Layer  Theory 
A  v  Experimental  Data  from  Fig.  17  in  Ref.  37 


Percent  X/C 

Fig.  19  Turbulent  Boundary- Layer  Parameters  on  a  Lifting  Infinite  Yawed  Airfoil 
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Fig.  21  Semilogarithmic  Plot  of  the  Turbulent  Boundary- Layer  Velocity 
Profiles  at  50-percent  Chord  on  a  Lifting  Infinite  Yawed  Airfoil 
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NACA  631-012  Airfoil  at  45.0-deg  Sweep 

Reoo,C  “  5,4  *  1q6*  Adiabatic  Wall,  =  0.74 

1  Present  Turbulent  Boundary-Layer  Theory 

q7DO  Experimental  Data  from  Fig.  16  in  Ref.  37 

Percent  X/C 

90  70  50  30 


Fig.  24  Streamline  Direction  across  the  Turbulent  Boundary  Layer  on  a 
Lifting  Infinite  Yawed  Wing 
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30°-Sweep  DAC  Airfoil,  C  -  L  0  ft 

Mqj  -  a  815,  pm  ■  3. 024  x  10*  Ibf/ft?,  ■  432°R 

Full-Scale  HIRT  Conditions 


X/C 


Fig.  25  Inviscid  Flow  Parameters  on  the  DAC  Airfoil  under  Full-Scale  HIRT  Conditions 
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Fig.  26  Adiabatic-Wall  Temperature  Distribution  on  the  DAC  Airfoil  under  Full-Scale  HIRT  Conditions 
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Three-Dimensional  Turbulent  Boundary  Layer 
30°-Sweep  DAC  Airfoil,  C-  1.0  ft 
MCO-0L815,  Re^c  "  L  05  x  108 
Full-Scale  HIRT  Conditions 
Adiabatic  Wall 

Constant  Tw  -  100a  0°R  (yT0  ^  -  2. 043) 
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Three-Dimensional  Turbulent  Boundary  Layer 

30°-Sweep  DAC  Airfoil,  C  =  1.0  ft 

Mcc  =  0.815,  Re„)C  =  1.05  x  108 

Full-Scale  HIRT  Conditions 

■Constant  T  =  1000. 0°R  (T  /T  =  2.043) 

Ay  W  O  f  oo 


x/c 


Fig  32  Stanton  Number  Distribution  on  the  DAC  Airfoil  under  ' 
Hot-Wall  Conditions 


fc' 
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Three-Dimensional  Turbulent  Boundary  Layer 
3CP-Sweep  DAC  Airfoil,  C  »  1.0ft 
Chordwise  Location  X/C  ■  0.90 


0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.0 


u  w 

ue  *  we 

Fig.  33  Hot-Wall  Effects  on  Turbulent  Boundary-Layer  Velocity  Profiles 
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Fig.  34  Hot-Wall  Effects  on  Turbulent  Boundary- Layer  Temperature  Profiles 
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TABLE  I 

NACA  63-1-012  AIRFOIL  ORDINATES  AND  INVISCID 
SURFACE  VELOCITY  DISTRIBUTION* 


~T/T, 

percent 

~77C, 

percent 

- 1 

:: 

Ui~ 

A  Ue,a 

V„  cos  A 

V.  cos  A 

Vqq  cos  a 

0.0 

0.0 

0.0 

0.0 

2.336 

0.50 

0.935 

0.750 

0.366 

1.695 

0.75 

1.194 

0.925 

r¥WWE  ?  wifi 

1.513 

1.25 

1.519 

1.005 

2.50 

2.102 

1.129 

0.933 

5.00 

2.025 

1.217 

fllPi  ^  *  ®gBi: 

7.50 

3.542 

1.261 

0.559 

10.00 

4.039 

1.294 

1.133 

0.484 

15.00 

4.799 

1.330 

1.153 

0.387 

20.00 

5.342 

1.349 

1.161 

25.00 

5.712 

1.362 

1.167 

30.00 

5.930 

1.370 

1.170 

0.249 

35.00 

6.000 

'1.356 

40.00 

5.920 

1.348 

1.161 

45.00 

5.704 

1.317 

1.148 

0.174 

50.00 

5.370 

1.276 

55.00 

4.935 

1.229 

1.109 

UMa  H  ft 

60.00 

4  .420 

1.181 

1.087 

0.121 

*5.00 

3.8*9 

1.131 

1 .063 

0.106 

70.00 

3.210 

1 .076 

1.037 

75.00 

2.556 

1.023 

1.011 

0.07O 

80.00 

1.902 

0.969 

0.984 

0.067  1 

85.00 

1.274 

0.920 

0.959 

0.055 

00.00 

0.707 

0.871 

0.933 

95. 0G 

0.250 

0.826 

C  .909 

100.00 

0.0 

0.791 

0.889 

0.0 

i 


Leading  Edge  Radius  =  1.037  percent  C 


*Table  taken  from  Abbott  and  von  Doenhoff  {Ref.  28,  p.  337) 
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APPENDIX  ill 

IMPLICIT  FINITE-DIFFERENCE  SOLUTION  OF 
GOVERNING  BOUNDARY-LAYER  EQUATIONS 


The  governing  turbulent  boundary- layer  equations  for  flow  over  an 
infinite  yawed  airfoil  as  derived  in  Section  2.  6  of  the  present  report 
form  a  set  of  coupled,  nonlinear,  parabolic,  partial  differential  equa¬ 
tions  with  the  transformed  coordinates  £  and  rj  as  the  independent  vari¬ 
ables.  For  completeness,  these  governing  equations  (Eqs.  (38),  (39), 
and  (40) )  are  repeated  below  with  primes  denoting  partial  differentiation 
with  respect  to  the  ^-coordinate: 

x  -  MOMENTUM 


£  f" 


GM 


f"  +  /9t0-{f')2] 


MOMENTUM 


(HI-  1) 


(III-  2) 


ENERGY 


fv(Q  *■  U«' 


+  Q'  = 


with  the  definitions 


{'  = 


U 


e 


(III- 3) 

(III- 4) 

(III- 5) 

(III- 6) 

(III- 7) 
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£  -  W 


£*  =  £[1  +-3 

V- 

r - .[>  .1 


fi 

SI 


dUe 

ue  ~df 


{?  (■  -  *;)  *  M)]pf  r'" 


(III- 8) 
(III- 9) 

(III- 10) 

(III- 11) 

(III- 12) 


The  boundary  conditions  are: 
MOMENTUM 


f(£r,  =  0)  =  0 
f'(£i?  =  0)  -  0 
lim  f  '(&jj)  =  1 

n~° 


(III- 13) 


c  (£lf  =  0)  =  0 

lim  c(f,Tj)  =  1 
TJ-wo 


ENERGY 


(III- 14) 


Hw  hw 
=  0)  =  —  =  — 

g  '(£*?  =  0)  =  0 

i1!  g(£>?)  =  i 


gw  (Prescribed  Wall  Enthalpy) 
(Adiabatic  Wall) 


(III- 15) 


Following  the  approach  by  Blottner  (Ref.  46)  and  Davis  (Ref.  47), 
the  momentum  and  energy  equations  (III-  1,  -2,  and  -3)  are  rewritten 
in  "standard"  form  for  a  parabolic  partial  differential  equation  as 
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d2W  dv  w  <?W 

i7+“'5?  +  a*f  *  ”457- 

0  (III- 16) 

where  W  =  V  for  the  x-momentum  equation,  W  =  c  for  the  z-momentum 
equation,  and  W  =  g  for  the  energy  equation.  Using  Eqs.  (Ill-  1,  -2, 
and  -3),  one  finds  the  coefficients  oi  through  a 4  in  linearized  form: 

x-MOMENTUM 

!fa? 

01  =  £* 

(III- 17) 

„  -<Sr 
“2=  f 

(III- 18) 

3 

(III- 19) 

-2£f' 

“4  -  g* 

(III- 20) 

z-MOMENTUM 

- 2 

(III-  21) 

a3  -  0 

(III-  22 ) 

a3  =  0 

(III- 23) 

a‘’  f 

(III- 24) 

ENERGY 

^(tt)  +  1  * 

- "  (9 

(III- 25) 

a2  =  0 

(III-  26) 
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(III- 27) 

(III- 28) 

The  rj  derivatives  in  Eq.  (Ill- 16)  are  replaced  with  finite- difference 
quotients  which  allow  variable  grid  spacing  in  the  17  direction  in  order  to 
concentrate  grid  points  in  the  region  near  the  body  surface  where  the  de¬ 
pendent  variables  change  most  rapidly.  The  derivative  in  the  ?  direction 
in  Eq.  (Ill- 16)  is  handled  in  the  usual  manner  as  a  two-point  backward 
difference  between  points  (m+1,  n)  and  (m,  n),  whereas  all  rj  derivatives 
are  evaluated  at  point  (im- 1,  n)  according  to  the  grid  mesh  shown  below. 


n' 


Station 


n  m+1 


The  solution  is  assumed  to  be  known  at  point  (m,  n)  and  unknown  at 
point  (m+1,  n)  so  that  the  finite-difference  scheme  to  be  constructed 
will  be  implicit  in  nature.  The  finite- difference  replacements  for  the 
derivatives  are  as  follows  (see  Appendix  IV  of  Ref.  13  for  derivation): 


*  *  V,  -  11  *  *>*J  m+-l 


(III- 29) 
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W]  IVi  - 


1 


rawl  Wm+l,D  -  ff- 

WELi.."  5? 


m  ,n 


(III- 30) 

(III- 31) 


where 


D1  =  (4n+i  ~  VD)  +  K2(j?n  -  i ?n_,) 

(Ill- 32) 

D2  =  <VH  -  +  K^n  "  %-\)2 

(III- 33) 

^n+l  ~  /  v 

K.  =  -  (constant) 

“  ^n—  1 

(III- 34) 

£  =  ij2  “  *?l  (constant) 

(III- 35) 

The  finite- difference  form  of  Eq.  (Ill- 16)  becomes,  upon 

substitution 

of  Eqs.  (III-29,  -30,  and  -31), 

AnWm+I.n+l  +  BnWm+l,n  +  CnWm+l,n-l  *=  Dn 

(III- 36) 

where 

-T  2  al 

(III- 37) 

"  °2  D1 

■B  -2(1  +  K)  a j(l  -  K2)  a4 

B°  =  D2  -  D,  +  °2  +  A? 

(III- 38) 

C  =  —  - 


2K  ^“l 


"  d2  U, 


a,  W_ 

4  m,n 


"ST 


(III- 39) 

(III- 40) 


For  Eq.  (Ill- 36)  to  be  linear,  the  coefficients  An,  Bn,  Cn,  and  Dn  must 
be  treated  as  known  quantities  at  point  n;  more  will  follow  on  this  sub¬ 
ject  later.  The  important  point  is  that  Eq.  (Ill- 36)  represents  a  set  of 
simultaneous  linear  algebraic  equations  under  this  restriction. 

Since  the  simultaneous  linear  algebraic  equations  resulting  from 
Eq.  (Ill- 36)  are  of  a  special  form  (tridiagonal),  an  efficient  method  of 
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solution  on  a  digital  computer  is  available  from  Richtmyer  and  Morton 
{Ref.  48,  pp.  198-201  and  274-282).  For  this  procedure  the  boundary 
condition  at  the  outer  edge  (n  =  N)  must  have  a  specified  value  Wm+i 
special  procedures  (see  Appendix  IV)  are  required  at  the  wall  (n  =  1). 
Because  of  the  special  form  of  Eq.  (Ill- 36),  the  relation 


exists  where 


W 


m-t-l.n 


E 


2 


e2 


E 


n 


e 


n 


E.  +  V  25»<N-1  (III-  41) 


“(A 2  +  C2  Fg) 


Bj+  Cj  Fj 

(III- 42) 

d2-c2  f, 

B2  C2  F2 

(III-  43) 

Bn  +  C„  En-1 

(III- 44) 

)  3  <  n  <  N-l 

a  ss 

“n  “  en— 1 

(III- 45) 

\  *  ~Cn  En-1 

with  the  numerical  constants  Fi,  F2,  and  F3  defined  as  in  Appendix  IV. 
The  quantities  E„  and  en  are  computed  from  Eqs.  (III-42)  through  (IH- 
45)  starting  with  n  =  2  and  progressing  to  n  =  N-l.  The  solution 
is  then  obtained  by  evaluating  Eq.  (Ill- 41)  from  n  =  N-l  to  n  =  2  with 
the  wall  value  (n  =  1)  evaluated  from  Eq.  (IV- 10)  of  Appendix  IV  accord¬ 
ing  to 


Wm+l,l  =  Fj  +  F2ffm+1,2  +  F3Wm+l,3 


(III- 46) 


where  F^,  F2,  and  F3  are  numerical  constants  defined  in  Appendix  IV. 

Knowing  the  distribution  of  the  f '  and  6  profiles  across  the  bound¬ 
ary  layer  from  the  above  procedure,  one  can  evaluate  the  transformed 
stream  function  f  from 


f<6?>  -/’ *'<£*>  *1  (III- 47) 

O 

where  the  integral  is  numerically  integrated  using  the  well-known  trap¬ 
ezoidal  rule,  viz. , 

=  .2,  [Hfrj  «  n)  +  nfrj -**_,)] D.  (Ill- 48) 
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with 


(III- 49) 


Inversion  from  the  transformed  (?,  57)  plane  to  the  physical  (x,  y)  plane 
is,  from  Eq.  (63), 


where 


(III- 50) 


(III- 51) 


because  of  the  constancy  of  static  pressure  across  the  boundary  layer. 
■Again,  using  the  trapezoidal  rule  method  of  numerical  integration  yields 

yn  =  2TU-  K  [6(^  "  ^  +  #6*  "  >71-1)1  Di  (HI'  52) 

r  e  u  e  1=  £ 

with  given  by  Eq.  (Ill- 49).  The  relationship  between  ?  and  x  is  deter¬ 
mined  from  integration  of  Eq.  (31) 


f-  (in- 53) 

O 

using  trapezoidal  rule  numerical  integration  to  advance  from  station  m 
to  station  m+1;  i.e.. 


Cn  -  (.  +  Kp *  (p.l-.D^,)^  (III- 54) 

where  Ax  is  the  x-direction  step  size  integration  increment  defined  by 

Ax  =  xm+ 1  -  xm  (III- 55) 

Note  that  arbitrary  x-direction  station  locations  are  allowed  via  the  a- 
bove  approach;  e.g,,  one  can  decrease  the  x-direction  step  size  in  re¬ 
gions  of  strong  adverse  pressure  gradient  in  order  to  accurately  define 
the  location  of  boundary-layer  separation. 

The  variable  n- direction  grid  mesh  used  in  the  present  work  is 
taken  from  Smith  and  Cebeci  (Ref.  49).  The  various  constants  used 
herein  are  as  follows: 

N  =  161,  K  =  1.063,  C  =  0.0010 
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Experience  with  varying  these  constants  and  observing  their  influence 
on  the  resultant  numerical  solution  has  indicated  that  the  above  choices 
are  adequate  under  the  present  subsonic  and  transonic  high  Reynolds 
number  flow  conditions.  These  values  may  not  be  satisfactory,  how¬ 
ever,  for  other  body  geometries  and  flow  conditions;  thus,  the  influ¬ 
ence  of  the  variable  grid  mesh  constants  should  be  ascertained  for 
each  new  investigation. 

The  mathematical  basis  of  the  above  tridiagonal  matrix  procedure 
applied  to  the  solution  of  boundary- layer  problems  is  from  Fliigge-Lotz 
and  Blottner  (Ref.  50).  The  present  application  differs  from  their  orig¬ 
inal  work  in  one  important  aspect:  the  linearized  difference  equations 
herein  are  uncoupled  and  solved  separately.  In  Fliigge-Lotz  and 
Blottner's  approach,  the  difference  equations  remain  coupled  and  re¬ 
quire  additional  machine  storage  and  manipulations  for  solution.  With 
the  present  uncoupled  approach,  the  difference  equations  are  iterated 
to  convergence  at  each  station  along  the  body;  thus,  one  must  pay  the 
price  of  iteration.  Through  this  procedure  the  linearizing  coefficients 
a\  through  (*4  are  evaluated  at  each  station  using  the  results  of  the  pre¬ 
vious  iteration.  Iteration  at  a  given  station  is  continued  until  successive 
values  of  f',  c,  and  g  differ  by  less  than  0.  10  percent  from  the  corre¬ 
sponding  value  of  the  preceding  cycle.  Typically,  about  three  to  eight 
iterations  per  station  are  required  for  both  laminar  and  turbulent  cases, 
with  the  number  of  iterations  per  station  increasing  as  the  pressure 
gradient  parameter,  j3,  increases.  Use  of  the  above  procedure  gives 
a  final  solution  at  each  station  which  is  exact  in  the  sense  that  it  repre¬ 
sents  a  converged  iterated  solution  to  the  governing  nonlinear  partial 
differential  equations  written  in  finite- difference  form. 

Along  the  stagnation  or  attachment  line  of  a  blunt-nosed  infinite 
extent  yawed  airfoil  where  ?  =  0  since  x  =  0,  the  governing  boundary- 
layer  equations  (III-l),  (III- 2),  and  (III- 3)  reduce  to  the  following  ordi¬ 
nary,  nonlinear,  differential  equations: 

x- MOMENTUM 


z- MOMENTUM 


pie  -  (f')2l 


s  0 


(III- 56) 


(III- 57) 
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ENERGY 


*:■ 

(&»•  • 


(W)«"  +  _^(?r)  +  fg'+  n' 


=  o 


(III- 58) 


with  the  boundary  conditions  (III- 13)  through  (III- 15).  In  order  to  ob¬ 
tain  starting  profiles  for  the  finite- difference  scheme  to  march  down¬ 
stream,  Eqs.  (Ill- 56),  (III- 57),  and  (III- 5 8)  are  solved  in  the  following 
manner  by  using  the  tridiagonal  matrix  procedure  described  previously. 
Along  the  stagnation  or  attachment  line,  ?  =  0  so  that  04  =  0  in  both  of 
the  momentum  equations,  as  well  as  in  the  energy  equation;  in  addition, 
the  term  containing  ?  in  ot\  of  both  the  momentum  and  energy  equations 
vanishes.  Furthermore,  jS  =  1  at  £  =  0  under  the  restriction  that 


near  x 


-•  -  a. 


=  0.  Initial  guesses  for  f',  c,  and  g  are  input  to  the  analysis  as 


f|'  =  1  -  exp  (-7) 


(III- 59) 


c 


I  = 


u 


(III- 60) 


gl  -  gw  +  U  -  8w)  f,'  (III-  61) 

where  the  subscript  I  denotes  the  initial  approximation.  The  equations 
are  then  iterated  to  convergence  in  the  same  manner  described  pre¬ 
viously;  an  averaging  scheme  is  used  to  speed  convergence.  Typically, 
about  10  to  20  iterations  are  required  to  generate  a  converged  initial 
solution  for  a  laminar  boundary  layer.  Approximately  60  iterations 
are  required  for  a  converged  turbulent  boundary-layer  solution. 


By  the  use  of  the  above  procedures  the  numerical  solution  of  any 
two-point  boundary- value  problem  governed  by  either  linear  or  non¬ 
linear  ordinary  differential  equations,  as  well  as  sets  of  coupled  para¬ 
bolic  partial  differential  equations  (either  linear  or  nonlinear)  is  re¬ 
duced  to  subroutine  status  on  a  digital  computer  in  that  only  the  coeffi¬ 
cients  q-i  through  04  must  be  defined  in  conjunction  with  the  required 
boundary  conditions  for  each  new  problem.  Based  on  experience  with 
analyses  of  this  type,  the  use  of  the  iterative  tridiagonal  matrix  approach 
where  applicable  is  highly  recommended. 
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APPENDIX  IV 

FINITE-DIFFERENCE  FORMALISM  OF  WALL  BOUNDARY  CONDITIONS 


Evaluation  of  heat-transfer  rate  and  shear  stress  at  the  body  sur¬ 
face  requires  numerical  determination  of  the  wall  derivatives  g'(£,  rj  = 
0),  c'(S,  rj  =  0),  and  f"(§,  r)  =  0)  since 

qw  “  g  —  °> 

-  '"«•!  -  « 

-  “lf.1-0) 


as  shown  in  Eqs.  (59),  (61),  and  (62).  Consistent  with  the  variable  re¬ 
direction  grid  mesh  spacing  defined  in  Appendix  III,  the  wall  derivatives 
are  evaluated  in  the  present  work  by  application  of  the  three-point 
Lagrangian  interpolation  formula  (Ref.  51,  pp.  71-77)  evaluated  ac¬ 
cording  to  the  following  diagram: 


(Point  n  = 


1) 
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which  yields,  at  station  m+1,  f 

f  "(,,=0)  =  £[B,(2  +  K)f"  +  B2(l  +  K)f'_  +  B-fLj  (IV- 1) 

c'(jj=0)  =  ^[B1(2  +  K)cJ?i  +  B2(l  +  K)c  +  Bgc,j^]  (IV-2) 

g'(»?=0)  =  +  +  B2(1  +  K)8,,2  +  (IV- 3) 


where 


and 


n  "I 

1  C2(l+K) 

(IV- 4) 

b2  =  +1 

2  KC2 

(IV-  5) 

B 

3  "  £2ku  +  k-) 

t 

(IV- 6) 

I'rj.  =  {'lv  = 

(IV- 7) 

O 

-3 

II 

O 

II 

-3 

(IV- 8) 

Br, .  *  6<7  -  l) 

(IV- 9) 

with  i  =  1,  2,  and  3.  The  parameters  K  and  £  are  defined  by  Eqs. 

(Ill- 34)  and  (III- 35),  respectively,  in  Appendix  III.  The  use  of  the 
three-point  Lagrangian  interpolation  formula  for  the  wall  derivatives 
is  consistent  with  the  formulation  of  the  variable  grid  spacing  deriva¬ 
tives  presented  in  Appendix  IV  of  Ref.  13  and  follows  the  work  of  Smith 
and  Cebecci  (Ref.  49,  Appendix  C). 


As  stated  in  Appendix  III,  the  tridiagonal  matrix  procedure  requires 
special  formalism  at  the  wall.  With  reference  to  the  above-presented 
sketch  (which  is  taken  to  apply  at  the  body  station  m+1  where  the  un¬ 
known  solution  is  to  be  determined),  assume  that  the  points  n  =  1,  2, 
and  3  are  related  according  to  the  general  relation 


Wj  =  Fj  +  F2W2  +  F3W3 


(IV- 10) 


where  Fi,  F2,  and  F3  are  numerical  constants  to  be  determined  and  W 
is  the  generalized  dependent  variable  defined  in  conjunction  with  Eq. 
(Ill- 16)  of  Appendix  III;  note  that 
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W1 

w2 

*3 


-  ^m+l,n=l 
=  ^m+l,n=2 


w 

m+ 1 ,  n=3 


(IV- 11) 


Using  Eq.  (TV- 10)  in  Eq.  (Ill- 36)  of  Appendix  III  evaluated  at  the  point 
n  =  2  yields 


+  1?2^2  +  C2(Fj  +  ^2^2  +  F3W3)  =  ^2 


(IV- 12) 


which  can  be  grouped  in  the  form 

W„ 


— (A2  +  CgFj) 


®2+  ^2^ 2 


P2  -  c2f, 
B2  +  ^2F2 


(IV- 13) 


Evaluating  Eq.  (Ill- 41)  of  Appendix  III  at  the  point  n  =  2  yields 


W2  =  E2W3  +  e2 


(IV- 14) 


which,  upon  comparison  of  Eqs.  (IV- 13)  and  (IV- 14),  results  in  the  re¬ 
lationships 


— (  a  2  +  n  2  Fj> 

B2+  C2  F2 


D2-C2F1 


e9  =  = - — 


®2  +  ^2  F2 


(IV- 15) 


(IV- 16) 


corresponding  to  Eqs.  (Ill- 42)  and  (III- 43)  of  Appendix  III. 


'  To  evaluate  the  numerical  constants  Fj,  F2,  and  F3  in  Eq.  (IV- 10), 
suppose  that  the  first  derivative  of  the  W  function  is  prescribed  to  be 
the  boundary  condition  at  the  wall  point  n  -  1.  Using  the  difference  ex¬ 
pression  given  by  Eq.  (IV- 1)  results  in 


+  B2n+K)R2  +  B3N3] 


Upon1  writing  Eq.  (IV-17)  in  the  form 


W 


1 


-B2(l+K) 
B  j(2+K) 


Wn 


Bj(2+K) 


W. 


CBj(2+K) 


(IV-17) 


(IV-  18) 


93 


AEDC-TR-73-1 1  2 


and  comparing  Eqs.  (IV- 10)  and  (IV- 18),  the  numerical  constants  Fj_, 
F2,  and  F3  become  defined  as 


W1 

-C(  1+K)W{ 

(IV- 19) 

j(2+K) 

2+  K 

-Bjd+K) 

(1+K)2 

(IV- 20) 

Bj(2+K) 

K(2+K) 

”B3 

-1 

Bj(2+K) 

K(2+K) 

(IV- 21) 

where  use  has  been  made  of  Eqs.  (IV-4),  (IV-5),  and  (IV-6). 


Now  suppose  that  the  value  of  the  W  function  is  prescribed  to  be 
the  boundary  condition  at  the  wall  point  n  =  1;  e.g.,  W-^  =  Ww.  For 
this  condition  Eq.  (IV- 10)  becomes 


Wi  -  Fi  +  F2W2  +  ^3^3  “  Ww 

which  requires  that 


(IV-22) 


(IV- 23) 
(IV-  24) 
(IV- 25) 


In  summary  of  the  above-derived  boundary  conditions  relative  to 
numerical  solution  of  the  boundary- layer  equations  via  the  tridiagonal 
matrix  procedure,  the  following  examples  illustrate  application  to  the 
wall  boundary  conditions  of  interest  in  the  current  work: 

x- MOMENTUM  EQUATION  BOUNDARY  CONDITION 
GIVEN  BY  EQ.  (Ill- 13) 


Set;  F  j  =  0 
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z- MOMENTUM  EQUATION  BOUNDARY  CONDITION 
GIVEN  BY  EQ.  (Ill- 14) 

Set:  F  j  =  0 


ENERGY  EQUATION  BOUNDARY  CONDITION  FOR  PRESCRIBED 
WALL  TEMPERATURE  gw  GIVEN  BY  EQ.  (Ill- 15) 

Set:  Fj  =  gw 


F3  =  0  -  I  : 

ENERGY  EQUATION  BOUNDARY  CONDITION  FOR  ADIABATIC  WALL 
GIVEN  BY  EQ.  (Ill- 15) 


Set:  Fj  =  0 

F2  Defined  by  Eq.  (IV-20) 

F3  Defined  by  Eq.  (IV-21) 

It  should  be  noted  that  the  above  formalism  allows  completely  general 
wall  boundary  conditions  which  can  vary  with  respect  to  location  on  the 
body  surface,  e.g.,  variable  wall  temperature  prescribed  over  the  for¬ 
ward  section  of  an  airfoil  with  the  aft  section  adiabatic. 
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